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PREFACE 


The Air Force SINCGARS VHF-FM radio is a frequency-hopping 
auci-jam radio which utilizes an electrically short antenna to 
minimize aerodynamic drag on airborne platforms. The development of 
optimally efficient, electronically tuneable antennas for this radio 
ts of interest. Although the antenna groundplane is 
platform-dependent, it is usually small compared to an rf wavelength. 
A circular groundplane provides a standardized groundplane geometry 
with which to model and evaluate candidate antennas. Accordingly, a 
VHF antenna range with an 8 ft. diameter circular groundplane has been 


constructed at The MITRE Corporation to evaluate candidate antennas. 


The electrical properties of a monopole element at the center of a 
circular groundplane of finite .adius are of interest to this program 
for (a) qualifying the antenna range; (b) establishing antenna 
standards with which to measure test antennas; and (c) modeling 
candidate antennas. Que survey of the literature revealed that 
although this antenna has the simplest geometry of any monopole 
antenna, its properties are neither well understood nor standardized, 
particularly for groundplane radii which are small or comparable to a 


wavelength. The-preeent paper attempts to address this deficiency. 
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SECTION 1 


INTRODUCTION 


Monopole antennas are commonly employed in airborne and ground-based 
communication systems at a wide range of frequencies. The electrical 
properties of such antennas are dependent upon the geometry of both the 
monopole element and the groundplane. Typically, the monopole element may 
be electrically-short (length much less than a quarterwave) or 
Near-resonant (length approximately a quarterwave) and may be thin (length 
to radius ratio much greater than 10°) or relatively thick (length to 
radius ratio of 10! to 10°). In addition, the groundplane dimensions may 
vary from a fraction of a wavelength to many wavelengths. Therefore it is 
desirable to know how the input impedance and radiation pattern of the 
antenna change as the dimensions of the monopole element and the 
groundplane vary. The directional gain on or near the radio horizon (the 
groundplane is assumed to be horizontal) is of particular interest since 
the maximum operational range of a communications system often depends upon 


the directivity on the radio horizon. 


This study is restricted to a monopole geometry consisting of a 
vertical cylindrical element at the center of a perfectly conducting, 
infinitely thin, circular groundplane in free space. This geometry is of 
interest because its radiation pattern is uniform in the azimuthal 
direction and because its electrical characteristics are a function of 
primarily only three parameters, namely, the element length, the element 
radius, and the groundplane radius, when each is normalized to the 
excitation wavelength. For these reasons this geometry is conducive to 


analysis, experimental verification, and standardization. 


A typical feed for the monopole antenna is a coaxial line whose inner 
conductor is connected through a hole in the groundplane to the vertical 
monopole element and whose outer conductor is connected by means of a 
flange to the groundplane. Typically, the inner conductor diameter is 
equal to the monopole element diameter and the outer conductor diameter is 
equal to the groundplane hole diameter. Unless stated otherwise, such a 
feed will be assumed in this study. The ratio of the coaxial line’s outer 
to inner conductor diameters affects the antenna’s input impedance, but only 
significantly for a relatively thick monopole element on a very small 


groundplane. 


For the idealized case of a groundplane of infinite extent and of 
infinite conductivity, the monopole antenna may be modeled by the method of 
images as a dipole with one-half the input impedance and double the peak 
directivity of the dipole. The infinite groundplane prevents monopole 
radiation into the hemisphere below the groundplane but allows a radiation 
pattern identical to that of the dipole in the upper hemisphere. However, 
for a monopole element mounted on a groundplane of finite extent, the outer 
edge of the groundplane diffracts incident radiation in all directions and 
consequently modifies the currents on the groundplane and vertical element 
from those of an infinite groundplane. At the outer edge of the 
groundplane, the currents on the top and bottom faces of the groundplane 
are equal in magnitude but opposite in direction because the net current 
must be zero at the edge. Outer edge diffraction becomes increasingly 
significant with decreasing size of the groundplane because of increasing 
magnitude of the currents on the groundplane faces at the outer edge. Edge 
diffraction can alter the input impedance by more than 100% and directive 
gain in the plane of the groundplane by more than 6 dB from the values for 


a groundplane of infinite extent. 


Theoretical models exist for predicting the effects of diffraction by 
the outer edge of the groundplane. The existing models may be classified 
into two categories distinguished by whether the current distribution on 


the monopole element is initially known or is unknown. 


When the monopole element is very thin and not too long, its current 
distribution is approximately sinusoidal and independent of the radius of 
the groundplane. Consequently, the element’s current distribution can be 
initially specified and only tke groundplane’s current distribution need be 
determined. For this category of monopoles, the theoretical models 
reported in the literature consist essentially of Bardeen’s integral 
equation method for the groundplane radius small compared to a 
wavelength‘), Richmond’s method of moments (groundplane only) for the 
groundplane radius not too large compared to a wavelength‘) | Leitner and 
Spence’s method of oblate spheroidal wave functions for the groundplane 


n(3I-(5) 


diffraction’© and the geometric theory of diffraction (GTD) for the 


radius comparable to a wavelengt Tang’s scalar theory of 
groundplane radius large compared to a wavelength, and Storer’s variational 
method for the groundplane radius very large compared to a 


(7),(8) 


wavelength 


When the monopole element is relatively thick, its current distribution 
is no longer sinusoidal and consequently the current distribution on both 
the monopole element and the groundplane need to be determined as a 
function of the groundplane radius. For this category of monopoles, the 
theoretical models reported in the literature consist essentially of 
Richmond’s method of moments for the groundplane radius not too large 
compared to a wavelength‘?) and Awadalla-Maclean’s method of moments 
(monopole element only) combined with the geometric theory of diffraction 
for the groundplane radius large or comparable to a wavelength(?) »(10) | 
Thiele and Newhouse have also reported a model which combines the method of 
moments with the geometric theory of diffraction but their computer program 


is unavatlable*'!?. 


Each of the existing models is valid for different and sometimes 
overlapping sets of limited values of groundplane radii. Some of the 
models require extensive numerical computation. For these reasons, the 
collection of models taken as a whole has several deficiencies. In the 
open literature, results for input impedance and directive gain have never 
been assembled as a continuous function of groundplane radius for the 
entire range of values from zero to infinite wavenumbers. In regions where 
models overlap, it is sometimes unclear which models are more accurate. In 
some models, numerical results have been reported for only a few values of 
groundplane radius. In one model (Bardeen’s integral equation) the base of 
the monopole element has not been allowed to be in the same plane as the 
groundplane which is the present case of interest. Computer programs are 
not available for some of the older models because the models were 
developed before the advent of computers. One of the most versatile of the 
models (Richmond’s method of moments) gives only the input impedance but 
not the radiation pattern. In one model (Leitner and Spence’s oblate 
spheroidal wave functions), one of the published algorithms for computing 
the eigenvalues is incorrect. Finally, extensive numerical results for 
small groundplanes and for resonant monopoles with finite groundplanes are 


not found in the open literature. 


This paper attempts to correct these deficiencies. Computer programs 
and numerical results are presented for all of the models. The induced emf 
method is utilized to determine the tnput impedance of a thin idealized 
monopole element in the absence of a groundplane. In Bardeen’s integral 
equation method, the excitation function for the groundplane currents is 
extended to include the singularity which occurs when the base of the 
monopole element is in the same plane as the groundplane. Richmond’s 
method of moments is extended to give the far~field radiation pattern. In 
Leitner and Spence’s method of oblate spheroidal wave functions, the 


published continued~fraction algorithm for computing the eigenvalues is 


corrected. Numerical results for input impedance and directive gain are 
presented as a continuous function of groundplane radius for arbitrary 
radius. Numerical results of various models are compared and the suspected 
best available results are identified. Extensive numerical results are 
given for small groundplanes and for resonant monopoles on finite 
groundplanes. New experimental data is presented and compared with 


numerical results. 


Circuit representations of the monopole antenna are developed in 
Section 2. Theoretical models and numerical results are presented in 
Section 3 for the case in which the current dlstribution on the monopole 
element is initially known. In Section 4, theoretical models and numerical 
results are presented for the case in which the current distribution on 
both the monopole element and groundplane are initially unknown. The 
theoretical models are compared with experimental data in Section 5. 
Computer printouts of directive gain and the far-field elevation pattern 
are given in Appendix A. Computer programs of the theoretical models are 


given in Appendix B. 


SECTION 2 


CIRCUIT REPRESENTATION 


2.1 Geometry and Coordinate Systems 


Consider a monopole element of length h and radius b which is located 
in free space at the center of an infinitely thin circular groundplane of 
radius a and of infinite conductivity (see figure 1). The groundplane 


radius, when expressed in wavenumbers, is given by 


(2.1-1) 


nm 
iit 
~ 
& 


where 

k = 27/ = wavenumber (m') 

A = excitation wavelength (m) 
The monopole element and groundplane have current distributions in real 
time given by 


I(z,t) = Re{1(zeI"], 1(p,t) = Re[1(pyeI”* 


] (2.1-2) 
where 
W = radian frequency of the excitation = 27c/) (rad/sec) 
c = wave velocity in free space = 2.9979 x 10° m/s 
I(z),I1(P) = element and groundplane current amplitude distributions, 


respectively (amp) 


A field point P(r, 8, 0), expressed in spherical coordinates, is shown 
in figure 1. The field is uniform in the azimuthal direction %. The 
relationships between spherical, cylindrical , and oblate spheroidal 
coordinate systems are shown in table 1. In the far-field, the elevation 
angle @ is related to the oblate spheriodal angular coordinate 7 by 


7 ~—~» cos 8 as the spheriodal radial coordinate & —+oo. 


P(r,6,) 


y 


Monopole Element at the Center of a Circular Groundplane 


TABLE 1 


Spherical, Cylindrical, and Oblate Spheroidal Coordinates 


SPHERICAL 
_ (2,0, 6) 


CYLINDRICAL 
(P, z, 0) 


OBLATE SPHEROIDAL 
(€,7,0) 


fatten : p= a{(1-n2)(14e2) 1/2 
elevation ) z=ané 
azimuth o 


Note 1: In Table 1, the notation 6,7 is that of Leitner and Spence (L+S), 
Franklin Institute Journal, Vol. 249, No. 4, pp. 299-321, April 1950. This 
notation is related to that of Abramowitz and Stegun (A+tS), “Handbook of 
Mathematical Functions”, National Bureau of Standards, Applied Mathematical 
Series, No. 55, p. 752, June 1964 by 


2 1/2 
its * Gays ~ 1) 


2 1/2 
Tas ~ CL ~Tass) 
Note 2: The cylindrical coordinates (p, z, %) are related to the rectangular 
coordinates (x, y, z) by 


x = P cos ) 
y= psin 6 
Z=2Z 


2.2 Directive Gain and Input Impedance 


At a far-field observation point P(r, 9, ), the numeric directive gain 
d(8,6) of the antenna is defined as the ratio of its radiated power density 
s(8,6) to its radiated power density averaged over all directions. The 
radiated power density averaged over all directions is equivalent to the 
power density radiated by a hypothetical “isotropic” antenna. Accordingly, 
the directive gain d(®,) expressed in spherical coordinates with the 


origin at the antenna, is given by 


s(8,0) 


27 7 
(1/4m s(8,0) sin 6 de db 
‘ee 


d(@,) = (2.2-3) 


For antenna patterns which are uniform in the azimuthal direction, such as for 


the antenna geometry of figure 1, Eq. (2.2-3) reduces to 


280) _ _2 s(n) 
fs) in @ de i -(n) an 
s sin s 
0 i 


d(8) 


= d°(7) (2 2-4) 


where 
2 
8(8) = (1/2) Jase, IHgl? = (1/2/eQ/ Ho IEg| 
Ry? Eo = far-field magnetic and electric field intensities, 
respectively 


s°(7), d°() = radiation power density and directive gain, respectively, 
in oblate spheroidal coordinates 


Jus € = wave impedance in free space 


The numeric directive gain d(@,6) is related to the directive gain D(@,) 
expressed in dBi by 


D(@,d) = 10 108, 4 d(9,0) (dBi) (2.2-5) 


The total time-averaged radiated power Pio of the antenna is given 


by 


tal 


10 


27 7 


Poi ie J f 8(@,0) r2 sine de dd 
0 “Oo 


7 
= 2m? fe eine de = 2m? [ s-(n) an (2.2-6) 
0 -1 

The antenna radiation resistance R, referred to the complex amplitude 
I(z=0) of the antenna base current, is defined by 

R=2P.__/|1(z=0)/? (2.2-7) 

total 

Substituting Eqs. (2.2-6) and (2.2-7) into Eq. (2.2-4), 


d(@) = s(@) 87r2/[R]1(z=0)|7] = 8°(7) 87r2/[RI1(z=0) 12] = d7(N) 
(2.2-8) 


The antenna input impedance, Z is given by 


in’ 


Z,. = V(z=0)/1(z=0) = Rit jx (2.2-9) 


where in 
v(z=0) = complex amplitude of the excitation voltage across the 
aperture of the coaxial line feed to the antenna (volts) 
R,. = input resistance (ohms) 
X,. = input reactance (ohms) 


The input resistance R is related to the radiation resistance R by 


in 
Rin =R+ Rate (2.2-10) 
where R is the ohmic loss resistance of the antenna for finite 


ohmic 
conductivity of either the monopole element or the groundplane. In the 


present paper, R = 0 because the monopole element and groundplane are 


ohmic 
assumed to be of infinite conductivity. Accordingly, 


Rin = R 


Equation (2,2-11) is a statement that the antenna is assumed to have an 


nT (2.2-11) 


efficiency of unity. 
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2.3 Relationship Between Radiation Resistance and Directive Gain 


on the Horizon 


For a vertical monopole element with a finite groundplane, the 
far-fiel' radiated power density on the radio horizon, s(9@ = 77/2), is 
determined only by the current distribution on the monopole element (and 
not the groundplane current) since only the vertical monopole element has a 
component of electron acceleration which is normal to the radio horizon. 
(This statement is not true for a groundplane of infinite extent since then 
a far-field point on the radio horizon is on the groundplane.) Identical 
monopole elements with tdentical current distributions, but mounted on 
groundplanes of different finite radii, will consequently have identical 
far-field radiated power densities on the radio horizon. Accordingly, 


3, (8 = 7/2) = § (8 = 7/2); identical monopole element, identical 
element current distribution, groundplane of finite extent 


where the subscript € denotes the radius in wavenumbers of the groundplane 
of arbitrary but finite extent (€ < 00) and the subscript 0 denotes a 


groundplane of zero extent (€ = 0). 


If Eq. (2.2-8) is substituted into Eq. (2.3-1) and the quantity 
s(6 = 7/2) |1(2=0) 12/8 re is computed, one obtains the following 
relationship between radiation resistance R and numeric directive gain on 


the horizon d(@= 77/2): 
d.(@ = 77/2) R, = d)(6= 7/2) Ry = constant; 


identical monopole elements, identical element current 


distributions, groundplane of finite extent (2,.3-2) 


12 


In Eqs. (2.3-1) and (2.3-2) the condition of identical element current 
distributions for groundplanes of different radii is generally not 
satisfied by monopole antennas. The element current distribution I(z) is 
generally dependent upon the groundplane current I(P) which in turn is a 
function of the groundplane radius. However, for monopole elements which 
are sufficiently thin electrically and not too long, the element current 
distribution I(z) is approximately sinusoidal and independent of the 
groundplane current I(P) (see Section 3.1). Expressions for d (= 7/2) and 
Ro» for elements with a sinusoidal current distribution, are determined in 
Section 3.2. Substitution of those expressions into Eq. (2.3-2), for the 
case of an infinitely thin monopole element (b=0), yields 


d(@ = 7/2) R = (N/AT) [1 - cos(kh)]7/sin@(kh); b=0, 
sinusoidal element current distribution, groundplanes of 
finite extent (2.3-3) 
where 
7 = wave impedance of free space = 376.73 ohms 


h= length of the monopole element 


The condition b=0 may be removed from Eq. (2.3-3) without substantially 
altering the result since the radiation pattern and radiation resistance of 
electrically thin elements, which are not too long, are weakly dependent 


upon the element radius (see Section 3.2). 
2.4 Characterization of Currents 


The characterization of the currents on the monopole element, 
groundplane, and coaxial line feed is illustrated in figure 2. The 
physical realization, circuit representation, and two circuit idealizations 
of the currents are shown in figures 2(a), 2(b), 2(c), and 2(d), 
respectively. In Fig. 2(b), the coaxial line excitation of Fig. 2(a) is 
replaced by equivalent electric and magnetic currents on a conductor 
completely enclosing the coaxial line. In Fig. 2(c), the normalized 
ferrlte attenuation distance hy/rA<< 1 is idealized to be zero. In Fig. 
2(d), the magnetic frill My is assumed to be negligible. 


13 


Lossy 


Figure 2. 


(c) 


(d) 


Characterization of Currents on Monopole Element, 
Groundplane and Coaxial Line Feed 


(a) Physical realization 

(b) Circuit representation 

(c) Idealization of circuit (hy,/A=0) 

(d) Idealization of circuit (h,/A=0, kb) « 1, 
I(z) or I(P) is initially known) 
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The currents of interest are the element current I(z) (positive in the 
+ z direction), the return current I (z) on the outside of the coaxial line 
outer conductor (positive in the + z direction), the current Tot on the 
bottom face of the groundplane (positive in the + P direction, and the 
current Top”? on the top face of the groundplane (positive in the 
negative P direction). A net groundplane current I(P) (positive in the 
positive $~ direction) is defined as 


TP) © Tige(P) = Tapp (P) (2.4-1) 

In the physical realization of the currents [figure 2(a)], lossy 
ferrite toriods are mounted along the outside of the coaxial line outer 
conductor. Such a procedure is commonly utilized on antenna ranges to 
reduce the radiation or pickup of unwanted signals from currents induced on 
the outside of the transmitter or receiver coaxial cables. For a 
sufficient number of ferrite toroids near the termination end of the 
coaxial line, the return current T(z) is approximately zero at distance 
greater than the l/e attenuation distance h, from the termination end of 
the line. In this paper, the ferrite toroids are assumed to be 


sufficiently lossy so that 


hy <r (2.4-2) 


where A is the excitation wavelength. 


The radii of the outer and inner conductors of the coaxial line are by 
and b, respectively, where b is also the radius of the monopole element. 
The wall thickness of the coaxial line outer conductor is assumed to be 
much less than its diameter. Consequently, the return current T(z) occurs 


at the radial cocrdinate ~P = b- 


The constraints on the various currents are: 


element I(z=h) = 0 (2.4-3) 
groundplane I(~=a) = 0 (2 .4-4) 
coaxial line I (2)* 0, -@O< z <h, (2.4-5) 
element-groundplane 1(z=0) = -I(P=b) (2 .4-6) 
groundplane-coaxlal line Tot P7>,) = I (220) (2.4-7) 


The element and groundplane constraints are a consequence of an open 
circuit at the end of the element and the groundplane. The coaxial line 
constraint is a consequence of the lossy ferrites. The element-groundplane 
constraint is a consequence of conservation of charge (Kirchhoff’s cuz~-ent 
law) at a node. The groundplane-coaxial line constraint is a consequence 
of conservation of charge along a conductor. 


(12) 


By the use of the equivalence principle the coaxial line feed 
excitation may be replaced by equivalent tangential field excitations 
defined along a surface completely enclosing the coaxial line. At field 
points external to this surface the equivalent field excitations will give 
the same field as the original source excitation. In the circuit 
representation of the monopole antenna currents [figure 2(b)], the coaxial 
aperture excitation is replaced by an equivalent surface magnetic current 
density (magnetic frill) My which sits on top of a thick groundplane of 
radius b> The magnetic frill My is defined to be positive in the positive 
azimuthal direction and is given by 


2 V(0)/{P1n(b,/b)], bSPSd, 
i : (2.4-8) 
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Eq. (2.4-8) is derived as follows. The radial field E of the coaxial 


line aperture, assuming a TEM mode excitation, is given by 


V(0)/P1n(b,/b), b< P< db, 
E, = (2.4-9) 
» PSb 


where V(0) is the positive voltage at z=0 across the aperture with the 
coaxial outer conductor at zero potential. By the equivalence 


(12) an aperture field may be replaced by a magnetic frill M 


principle 
which sits on top of a groundplane congruent with the aperture surface and 
which is defined as 


— — 


(2.4-10) 


M = 


xn 
tangential 


where n is the outward normal to the aperture surface and E. is 
tangential 
the tangential field at the aperture surface. Substituting Eq. (2.4-9) 


into Eq. (2.4-10), 
= nN nN n~ n~ 
Me (uy E,) xuls- Eo = ug (2.4-11) 


where Mg is given by Eq. (2.4-8). Eq. (2.4-8) agrees with the result 


obtained by Richmond 2), 


In the circuit representation of figure 2(b), the net groundplane 
current I(P) is the same as defined by Eq. (2.4-1) with the additional 


current constraint 


Tort?) 70, OPK d, (2.4-12) 
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If the circuit representation of the monopole antenna in figure 2 is 


now idealized by setting the ferrite 1/e attenuation distance h,/A=9, then 
I (2) =0 for z <Q, h,/ A= 0 (2.4-13) 


Consequently, the coaxial line outer conductor may be removed from the 
circuit as shown in figure 2(c). The groundplane — coaxial line current 
constraint of Eq. (2.3-7) is not disturbed by such an idealization. Since 
it has already been assumed that hi A, the idealization h,= 0 does not 
appreciably alter either the radiation pattern or the input impedance of 
the monopole antenna provided that the monopole length h >> h,- Finite 
currents, on an aperture which is small compared to the excitation 
wavelength, contribute little to the far-field and input impedance 
(compare with the results for an electrically-small dipole). Experimental 
radiation patterns and measurements of input impedance (see Section 5) 
confirm that the use of lossy ferrite toroids along the coaxial line outer 
conductor yields results which are in close agreement with theoretical 
results for the circuit idealization condition of Eq. (2.4-13) -- even for 


electrically~small groundplanes. 


In the idealization of the monopole antenna circuit, the magnetic frill 
My may be removed [figure 2(d)] without appreciably altering the radiation 


pattern or input impedance, provided that 


kb, << 1 
conditions for neglecting magnetic frill 
I(z) or I(p) is initially known (2.4-14) 
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The condition kb, << 1 corresponds to the condition for TEM mode excitation 
of the coaxial line and for negligible power being radiated from the 
coaxial line aperture. If either I(z) or I(P) is initially known then the 
one which is not known may be determined from the other without requiring a 
knowledge of the original coaxial line excitation or its equivalent ~- 
provided that the field radiated by the known current distribution is the 
predominant field incident on the conductor whose current distribution is 
unknown. When neither I(z) nor I(P) is known, then the original source 
excitation or its equivalent (in this case, the magnetic frill My) must be 
specified to determined the unknown radiated fields. 
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SECTION 3 


MODELS IN WHICH THE CURRENT DISTRIBUTION ON THE MONOPOLE 
ELEMENT IS INITIALLY KNOWN 


3.1 Boundary Conditions 


The current amplitudes, I(z) and I(P), on the monopole element and 
groundplane, respectively, are generally complex and initially unknown 
quantities (see figure 1). Consider now the case where the current 
distribution on the monopole element is assumed to be sinusoidal. For such 


a case and for the waveform dependence given in Eq. (2.1-2), 


1(0) 


Mz) = Stach) 


sln{k(h-z)], O¢ z¢h (3 .1-1) 
where 
h = monopole length (m) 
k = 27/X= wavenumber (mt) 
1(0) = current amplitude of the monopole base current at z=0 (amp) 
From Eq. (3.1-1) 
arg I(z) = constant, O¢ z¢h (3.1-2) 
Although a sinusoidal distribution of current is not possible even for 


an infinitely thin antenna, Eq. (3.1-1) is most likely a fair approximation 
to the current if the monopole element is sufficiently thin electrically 
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(13) 
(14) 


and not too long. 


length 2h, Elliot 


For a center-fed dipole of radius b and total 
gives examples where the current distribution is 
approximately sinusoidal and is of approximately constant argument for 
b/A= 1.0 x 10” and h/ = 0.125, 0.25. However for b/A = 1.0 x 104 and 
h/A = 0.375, 0.5 Elliot demonstrates that the current distribution is no 
longer sinusoidal near the center of the dipole nor is arg z approximately 
constant. Balanis‘!5) shows that for b/A= 2.0 x 107" and h/A = 0.25, 0.5 
the current distribution is not sinusoidal near the center of the dipole. 
2 and h/A = 
0.125, 0.25, 0.375, 0.50 the current distribution is neither sinusoidal nor 


Elliot and Balanis also demonstrate that for b/A = 1.0 x10 


of constant phase and that the deviations from Eqs. (3.1-1) and (3.1-2) 
increase with increasing values of h/A and b/X. On the basis of the above 
results, it appears that Eqs. (3.1-1) and (3.1-2) are approximately valid 


for the conditions 
b/A< 10%, h/AK 0.25 (3.1-3) 


In addition to the constraint on I(z) given by Eq. (3.1-1), assume that 


the return current I (z) on the outside of the coaxial line outer conductor 


(see figure 2) is given by 
1(z)®0 ,2z<¢ 0 (3.1-4) 


The constraint of Eq. (3.1-4) corresponds to the idealized condition that 


the ferrite toroids have a 1/e current attenuation distance h given by 


hy/A= ) (3.1-5) 


It should be noted that Eqs. (3.1-4) and (3.1-5) do not alter the 
constraint I (220) = Thor O7 b,) given by Eq. (2.3-7) nor do they impose 
any constraints on Tot 2 b)) where Thor P* b)) is the current on the 


bottom of the groundplane at a radius equal to that of the outer conductor. 
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Combining the current constraints given by Eqs. (3.1-1) and (3.1-4) 
with those given by Eqs. (2.4-3) - (2.4-7) and Eq. (2.4-12), the current 


constraints on the monopole antenna are given by 


element I(z) = [1(0)/sin(kh)]sin[k(h-z)], 0 £z¢h (3.1-6) 
groundplane I(@=a) = 0 (3.1-7) 

Thor (P)? = 9% O<PK< db, (3.1-8) 
coaxial line I (2)=0, z<o (3.1-9) 
element-groundplane I(z=0) = -I(=b) (3.1~-10) 
groundplane-coaxial line Tot (?) = I (z=0) (3.1-11) 


In Section 3, it will be assumed that all models satisfy the current 
constraints given by Eqs. (3.1-6) - (3.1-11). The results are expected to 
be approximately correct if the monopole element is electrically 
sufficiently thin and not too long [Conditions (3.1-3)] and if the ferrite 
toroids are sufficiently lossy [h)/A<< 1, idealized by Condition (3.1-5)]. 
For these conditions, the circuit representation of the monopole antenna is 
shown in figure 2(d). 


For the current constraints of Eqs. (3.1-6)-(3.1-11), the total 


magnetic and electric field intensities H¢total) —¢(total) 


at an arbitrary 
field point P(x,y,z) external to the element excitation source points are 
simply the vector sum of the fields resulting from the element current and 
the current induced on the groundplane by the fields incident from the 


element. 
Accordingly, 


current constraints of 
» Eqs. (3.1-6)-(3.1-11) (3.1-12) 
(e) Ho} ? 


“ltotal) _s(e) 4 8) FEE) ACE) 


~p(total) _-E(e) + 8) TE 
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where 

He) Ele) = magnetic, electric field intensities, respectively, generated 
by the element current 

“6 E62) = magnetic, electric field intensities, respectively, generated 
by the groundplane current induced by the element incident 


fields. 


The element fields He”, EO 


groundplane fields H¢® , “E¢2) 


are determined in Section 3.2. The 
are determined in Sections 3.3 - 3.8 for 


groundplane radii of various extents. 
3.2 Induced EMF Method, Groundplane of Zero Extent 
Concept of a groundplane of zero extent 


Consider a monopole antenna excited by a coaxial line whose outer 
conductor of radius by is terminated by free space rather than by a 
groundplane [figure 3({a)]. The groundplane for such an antenna is denoted 
as being of “zero extent”. As was shown in Section 2.4, the coaxial line 
excitation may be replaced by an equivalent magnetic current (frill) My 
sitting on top of a thick groundplane of radius P= by {figure 3(b)]. For 
sufficiently lossy ferrlte toroids along the outside of the coaxial line, 
the current on the exterior of the coaxial line outer conductor may be 
neglected [figure 3(c)]. The magnetic frill may be removed from the 
circuit without appreciably affecting the results since kb, << 1 for the 
assumed sinusoidal current distribution on the monopole element [see Eq. 
(3.1-3) and the discussion concerning the circuit idealization of figure 
2(d)]. Finally, the groundplane of radius by may be removed from the 
circuit without appreciably affecting the results since a finite current on 
an electrically-small conductor does not radiate appreciable power compared 


to the power radiated by a monopole element of length h >> b The circuit 


1° 
idealization of a monopole antenna with a groundplane of zero extent is 
therefore an electrically thin monopole element with no groundplane [figure 


3(d)]. 
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Characterization of Currents for Monopole Antenna with Groundplane of Zero Extent 
(a) Physical realization 


(b) Circuit representation 
(c) Idealization of circuit (h,/A=0) 
(d) Idealization of circuit (hy/A=0, kby «1, I(z) or I(P) is initially known) 


The near-fields, far-fields, and input impedance of an electrically 
thin monopole element are derived and summarized in the remainder of 


Section 3.2. The input impedance is derived by the induced emf method. 


Near-fields 


Consider a monopole element of length h with a sinusoidal current 
distribution I(z“) = [1(0)/sin(kh)] sin [k(h-z~)]. 0 < 2° ¢h, at points 
Q(x”, y”, 2°) on the surface of the element (figure 4). For an 
electrically-thin element with a known current distribution, the fields at 
arbitrary points P(x, y, z) external to the element may be determined 
almost exactly be approximating the source points to lie on the element 
axis, i.e., Q(x”, y~, 2°) ~Q(0, 0, z”). For the current waveform of 


Eq. (2.1-2), the magnetic vector potential A is given by‘) 


2 1(0) h _ 
A(x,y,z) = vA, = a, tay Jf sinjk ae exp(-ikP eer 
where ~7 
uy” ee eae of a space = 47x10 henries/m 
BQ = [(x-n")? + (y-y")? + (ae alg 
ise z B ules yy _ (p? ri (2-27)? 1? 


nN NOUA 
Up» Yoru, = unit vectors in the radial, azimuthal, and axial 


cylindrical directions, respectively, 


The magnetic and electric fleld tntensities, H and E, respectively, are 


given by 


- = 1 (9 A) = ies 1 _OAg 
Ho . 6 Ho ~ yp (3.2-2) 
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P(x, y, Z) 


(p, 7712,4) 


Figure 4. Monopole Element Geometry 
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= 1 = | Jy 1 
jue, VED" Np joe “Se * Ye jue 3p 


-12 
where &. * permittivity of free space = 8.854 x 10 farad/m 


—_ — —-~z 
Exact closed form expressions of A, H, E and radiation resistance were 


7 
first obtained for an infinitely thin element by Brillouin‘! ) 


(18) 


and are 


summarized by Stratton The magnetic and electric fields are given by 


=J2 sin(kh) ex(-3ee)| (3.2-4) 
p= Lorn | C0) expC Skee) ig cosh) exp(-Jkr) 
p  —- 40sin(kh) r . 


y sin(kh) a | te exa(-see)|| 
Zz r 
j.1(0) 7 (2rh) “expC Ike.) _ 2  cos(kh) exp(- jkr) 
4x0 sin(kh) vO r 
2 2 
, sin(kh) sere) oats #*]| (3.2-5) 
kr r L 
ER = -j 1(0)7 exp(—Jjkr,) _ cos(kh) exp(- jkr) 
z 47 sin(kh) tS r 


4 sin(kh) cs exp(- jkr) 
Oz r 


_ _-4 10) | exp(“Jkro) cos(kh) _exp(-Jkr) 


47 sin(kh) tO r 


_ 2 sin(kh) er ee + = ]| (3.2-6) 


r kr 
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where 
= k/We = wave impedance of free space = 376.73 ohms 
22 2-1/2 
rot (Po + (2-h)*] 
ew (p2 Fs 22) 1/2 
Eqs. (3.2-4) - (3.2-6) are identical to the results given by Stratton 
(& =h, 2 #0, 2=h, Oskhy rye ty ey ey 1s -1(0)/{sin (kh)] 
after the substitution of -j for j to account for the exp (-jwt) waveform 
of Stratton instead of the exp( jut) waveform of Eq. (2.1-2). The fields 
given by Eqs. (3.2-4) - (3.2-6) are exact for an infinitely-thin element 
and are almost exact for an electrically-thin element with the same 
sinusoidal current distribution. 
Far-fields 
Consider a field point P at a sufficiently large radial distance r 


which satisfies both the far radiation zone and Fraunhofer diffraction 


conditions given by 


h << r, kr >> 1; far radiation zone conditions (3.2-7) 
kh? /2r <€ 27; Fraunhofer diffraction condition (3.2-8) 


For these conditions, 
exp(-Jkr,)a exp[-kj(r-h cos®)], h << r, (kh?/2r) << 27 
l/r j= l/r, hr 
z-hsxz = rcos®, hh «<r (3.2-9) 
(l/kr) - (22/kr3) - (j22/r2) x ~-3 cos*®, kr >> 1 


(4/r?) + (1/kr3) x ( 4/4"), kr >> 1 
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For the approximations of Eq. (3.2-9), the “exact” fields given by 
Fqs. (3.2-4) - (3.2-6) reduce to the far-fields given hy 


_ _j 1(0) exp(-ikr) t P a _ 
ar = iar sin® sin(kh) exp (jkhcos@) cos(kh) j cosA sin(kh) 


(3.9-19) 
j_1(9)cos® exp(-jkr) [ | 
Eo eS OOS |1:SO-s«Scos(kh) - { cos®@ sin(kh) 
fe) 47rr sin® sin(kh) (3.9-11) 
= 0 = 
E = oe ~ cos(kh) - j cos@ sinh) 


(3.9-17) 
aa nw 
The resultant electric field F =unF, + uF. reduces in the far field to 


E = Oye g> far-field (32-13) 


where 


4 100) exp(-ikr) 


Fa = —Amtr sin® sin(kh) [; ad cos(kh) ae j cos9 sinckh) 


—= —_>_ om 
The time-averaged Poynting vector S = (1/2) (F x H) = 


. nN 
(1/2) (uy Eg x ug Hy = a s(9) with a time-averaged radiated power densitv 
s(9) given hy 


> 
s(@) = IRQ !°/(2) =(0/7) IHgl” = EO (3.9-14) 
3°97 r sin’ (kh) 
where 
2 9 
erqy= leosGn cos8) - coy + [sin(kh cos@) - cos 9 sin(kh)] 
sin 9 


The direction of maximum radiated power is @ = 7/2 rad. 
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Substituting Eq. (3.2-14) into Eqs. (2.2-6) and (2.2-7), the radiation 
resistance R, referred to the base-current, is given for an infinitely thin 


element by 


is 
7 f £(98) sin® de 


7 (ohms ) 
87 sin’ (kh) 


m 7 [Cin(2kh) - sin@(kh)}, b=0 (3.2-15) 


47M%sin (kh) 


dt 


Zz 
where Cin(z) modified cosine integral ©)? ). f = (l-cos t). 
0 


The result given for the definite integral of Eq. (3.2-15) is readily 
obtained by letting t = cos 9 and noting that 


1 1 [ 1 | 
= — — >- ——— 
1-2 2 x y 


where x = 1l+t and y = l-t. 
Substituting Eqs. (3.2-14) and (3.2-15) into Eq. (2.2-8), the numeric 


directive gain d(@) is given by 


£(9) 


d(@) = 7) 
Cin(2kh) - sin (kh) 


(3.2-16) 


where £(9) is defined in Eq. (3.2-14). 
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For electrically-short monopole elements (kh << 1) and for = 376.73 
= 1207 ohms, Eqs. (3.2-14) - (3.2-16) reduce to 


4 2 
£(@) = + (kh) sin® an << 1 (3.2-17) 


15 11(0)17 (kh)? sin2e 


s(@) = 5 , kh (3.2-18) 
16 Wr 
R = 5 (kh)? (ohms), kh << 1, b = 0 (3.2-19) 
3 si 29 
d(@) = nv kh <1 (3.2-20) 


The relative power radiation pattern 8(@)/s = s(@)/s(m/2), radiation 
resistance R, and the directive gain (directivity) d(@) are tabulated in 
Appendix Al, Tables Al-1, Al-14, Al-26, Al-38 for h/A = 1/4, 1/10, 1/40, 
and 1/100, respectively. The numeric directive gain pattern d(@) for a 
quarterwave monopole element [see figure 8(a) in Section 3.9] is similar to 
that of a half-wave dipole except that its peak directivity is less (1.88 
dbi vs. 1.76 dBi) and its 3 dB beamwidth is more (94 deg vs. 78 deg). 


Input Impedance 


The input impedance may be determined by the induced emf method 


introduced by Brillouin‘!”), With reference to figure 4, the input 
impedance rr is given py29) 
h 
Zz ue = = a = = “a a 
in at) fi 1,07 b, z 2°) ECP b, z= 2°) dz”, b> 0 
0 


(3.2-21) 


where 


a = 1(0)/sin(kh). 
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The induced emf method is indeterminate for collinear current elements 
unless the elements are of radius b > 0. For sufficiently thin elements, 
the current L is given by Eq. (3.1-1) and the electric near field E. is 
given almost exactly by Eq. (3.2-6). Substituting Eqs. (3.1-1) and (3.2-6) 
into Eq. (3.2-21), 


exp(-jkr.) 


h 
= Ui , =27)}{-——___ —_* 
Zig” rarateny Jf 4° statkdnrs | r, 
0 


_ cos(kh) _exp(-jkr) a sin(kh) 3 _ [ exp(-jkr) 
r Oz* r 


J -20 


where 
[b24(27-h) 21/2 


(b2 oe 2-2yt/2 


ta] 
] 


ta] 
W 


Each of the three terms of the integrand of Eq. (3.2-22) may be integrated 
by the methods summarized by Stratton‘!®), However, the third term of the 
integrand should first be integrated by parts. 


Accordingly, 


Zan = Rye + j Xin = Z, + 25 + Z3 (3.2-23) 
where 
h 
z* LJ J dz” sin[k(h-z)} S*P(-dkr,) 2 R, + 5X, 
4xein? (kh) r 
0 o (3.2-24) 
R, = - “5 : Cin(x,) + + Cin(x,) - Cin(x,) 
4nsin (kh) 
(3.2-25) 
7) 1 1 
4xsin’™ (kh) 
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2 j ul -dz° sin[k(h-z~)] cos(kh) exp(~ jkr) = Ro+ 5X, 
4nsin’ (kh) 0 


(3.2-27) 
n - sin(2kh) 
= ———_ | —— _ [si(x,) - Si(x,)] 
Angin? (kh) : ? 
fede’ (hd (CIN eS Cla) FS ce | (3. 2-28) 
3 2 2 2 1 
= n - gin(2kh) e 
eg, | EE timex cin(x,) + £n(x,/x,)] 
2 (kh) 
4: £08 eh) tsi (x,) + Si(x,) - 2 six) (3.2-29) 
h 
= it fe sin{k(h-z')] sin(kh) 5 [eset | R,t+ 5X, 
4tsin (kh 0 (3.2-30) 
2 
= ——- ~£40 EN) Leena, ) + Cin(x,) ~ 2 Cin(x,) 
4ngin‘ (kh) 2 
: 2(x, + x,) (sin x) + sin x) : 4 X4 sin x4 
x + x5 + 2(kb)* x5 + (kb)? 
2(x, + x,) (cos x, - cos x,) 
. stn(2ih) $i(x5) z Si(x,) Z 1 : 2 5 1 , 2 
x] + Xo + 2(kb) 
(3.2-31) 
2 
n sin (kh) 
= ——__— | St" Isi(x,) + Si(x,) - 281(x,) 
- 2c) a i 4 a 
2(x, + Xo) (cos x) + cos Xy) 4 X34 COS X4 ! 
+ A a RE NS RE RR OIE a NC a ee cE a 
x + x5 + 2(kb)” x5 + (kb) 


+ —Sint2kh) | ~Cin(x, ) + Cin(x,) + Ln(x,/x,) 


2(x, + Xy) (sin x, 7 sin xy) 
Oty ot ts ee gta ae 2e) 
x + Xp + 2(kb) 
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where 


x > kL (b + h2y!/2 4 ny 
2.1/2 
x, ® k[(b? + h2y)/2 - hy 
X4 = kb 
x 
Cin(x) = modified cosine integral‘??? -f a (l-cos t) 
0 
(19) * sine 

Si(x) = sine integral = I gee dt 


Summing Eqs. (3.2-23) - (3.2-32), 


feta, + Cin(x,) - 2 Cin(x,) 


R,. = ——>— 
in 4nsin? (kh) 
P sin(2kh) E + x, ) (cos x, 7 cos x5) 
2 2 2 ri 2 
x) + Xo 2x3 
2 (x, x») (sin x) + sin Xo) sin X3 
+ sin’ (kh) | ———>——__—_____—_ - ——_ 
2 2 2 x 
x) + Xo + 2x, 3 
(3.2-33) 
XxX, = cae eee Si(x,) + Si(x,) - 2 Si(x,) 
in 4xgin® (kh) 1 2 3 
+ - 
__sin(2kh) (x) x, ) (sin x) sin x,) 
2 2 
3 


2 2 
xy + Xo + 2x 


+ sin’ (kh) je i. x 
1 2 


(x, + x, )(cos X, + cos X») cos xX, | 
x +x + 2x 


(3.2-34) 
where x, = k{(b? + n2y!/? + hy, x, = kt(b? + n2y!/? - ny, x = kb. 
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For b=0, Eqs. (3.2-33) and (3.2-34) reduce to 


R, = —_t Cin(2kh) - sin-(kh) » b=0 (3.2-35) 
" 4nsin* (kh) 
n sin? (kh) 
Xi. a __ | Sa (2khy) - SPY |, bo (3.2-36) 
4nsin (kh) 


Eq. (3.2-35) agrees with the result for radiation resistance given by 
Eq. (3.2-15). A comparison of Eq. (3.2-33) with Eq. (3.2-35) reveals that 
the input resistance is relatively insensitive to the monopole element 
radius b for kb << 1. The input resistance given by Eq. (3.2-35) is 
plotted in figure 5. For kh =7/2 radians, Ree = 19.4 ohms. 


The input reactance is sensitive to the element radius b as seen in 


Eq. (3.2-36). For b = 0, the input reactance is given by 


- co ohms, kh # n7, n=1, 2,3... 


x ,b=0  (3.2-37) 


ooohme, kh =n7m, n=l, 2,3... 


From Eq. (3.2-36), resonance (KX = 0) occurs for 


(kh) = n+ [kb si(2nm})/?, 


n=1, 2,3... (3.2-38) 
resonance 


The input resistance at resonance, for kh given by Eq. (3.2-38), is found 
from Eq. (3.2-35) to be 


oF n cin(2n7) 


ORG mabaanee ~ 47 kb Si(2n7) 


- 1 (3.2-39) 


Minima of input resistance occur for kh given by 


(kh) Oe ee ene | es ne ee ee (3.2-40) 


minima 
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RADIATION RESISTANCE, Ry (OHMS) 


1) =376.73 OHMS 
si ¥=0.57771.... 


ie) 0.577 1.077 1.577 2.077 2.577 3.0 °T 3.577 


NORMALIZED MONOPOLE LENGTH, 2m h/A (WAVE NUMBERS) 


Figure 5. Radiation Resistance of Electrically-Thin (kb]1) 
Monopole Elements 
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The input resistance at these minima element lengths is found from 


Eq. (3.2-35) to be 


oes 
minima — 47 


(R,_) [y- 1+ £n[(2N + 1)7) (3.2-41) 


Bi 
where 


Y= Euler’s constant = 0.57721 


The accuracy of Eq. (3.2-40) increases as N —*> ©. 
Summary of Results 


The input impedance and directive gain properties of quarterwave and 
electrically-short monopole elements with groundplanes of zero, large, and 
infinite extent are compared in Table 4 of Section 3.9. The peak 
directivity is approximately 3 dB less with groundplanes of zero extent 
than with groundplanes of large extent. However, the directive gain on the 
horizon is approximately 3 dB more with groundplanes of zero extent than 
with groundplanes of large but finite extent. The radiation resistance 
with groundplanes of zero extent is approximately one-half that with 


groundplanes of large extent. 


Unlike dipole elements in free space whose first resonance occurs for 
dipole half-lengths approximately equal to a quarterwave, monopole elements 
with groundplanes of zero extent have a first resonance for an element 


length approximately equal to a half-wave. 
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3.3 Integral Equation, 0 <¢ ka ¢ 2.75 


In Section 3.2 the fields generated by the monopole element were 
determined. These fields impinge on the groundplane and induce a 
groundplane current. For sufficiently small groundplanes, the fields 
generated by the groundplane current may be determined by Bardeen’s 


integral equation method"? 


In this method a cylindrically symmetrical electromagnetic wave 
(generated by the element) is incident on the groundplane disc. The 
fields generated by the induced groundplane current are required. The 
solution depends upon solving an integral equation of the first kind. 
For groundplane radii of arbitrary radius, the integral equation is 
not readily solvable because it contains two integrals. However, for 
sufficiently small groundplane radii, Bardeen neglects one of the 


integrals so that the integral equation may be solved explicitly. 


Although Bardeen gives a general formulation of the solution for 
the resulting single tntegrand integral equation,his only explicit 
results are for the case when the incident wave is generated by an 
infinitely thin dipole element whose base is at a non-zero height 
above the center of the groundplane. Bardeen restricted his solution 
to elements at a non-zero height above the groundplane {in order to 
avoid having a source point (the base of the element) at a near-field 
point of interest which, for the integral equation method, includes 
the center of the groundplane. The total field is then given by the 
vector sum of the incident field and the {induced field [see Eq. 
3.1-12]. 
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The present case of interest is a monopole element whose base is 
in the plane of the groundplane, f.e., at a vertical height v=0 above 
the center of the groundplane (see figure 6). For this case, the 
total magnetic field may be determined by first evaluating the ffeld 
with the element at an arbitrary height v > 0 and then by evaluating 
the resulting expression in the lim v-0. By such a procedure, an 
indeterminate expression is avoided for the field generated hy the 


induced groundplane current. 


(total) 


Accordingly, the total magnetic field intensity ar [see Fq. 
(3.1-1.2)] in the limit v-—0 is given by 
vim H6EOCAL) yam Hoe) 4 1am Af? (3.31) 
v—0 v—0 v0 
(e) (g) 
where Ho and Hg are the magnetic field intensities generated by 


the element and groundplane currents, respectively. In the following 
evaluation of the two terms of Fq. (3.3-1), the fields are assumed to 
have an elt time dependence [see Fq. 2.1-2] unlike the en dte time 


dependence assumed by Bardeen. 


The first term of Eq. (3.3-1) is given exactly, for an infinitelv 
thin element, by [see Fq. 3.2-4] 


(e) . j 1(0) _ _ - 
oo Uap sin(kh) [exo jkr) cos(kh) exp(-ikr) 


iz sin(kh) exo(-th)| 


(3.3-7) 
r 
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P(p,d,z)= P(r 6,6) 


= 


.) 
a 
\\ 


Figure 6. Thin Monopole Element Whose Base [Is Above the Center 
of a Circular Groundplane 
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In the far-field, Fq. (3.3-2) reduces to [see Fq. (3.2-10)] 


lim ny”? = A [4 exp(ikh cos®) + 4 cos(kh) + cos® sin(kh)] , 
v—0 ? far-field (3.3-3) 
where 
eo & Ge 
ro) 4mr sin8 sin(kh) ~ 


The second term of Eq. (3.3-1) may be evaluated by utilizing 


(2) 


Bardeen’s Eq. (31) for ur) which he obtained as a solution to his 


single integrand integral equation. For an infinitely thin element 


and for sufficiently small groundplanes (ka<1), ae is given by 


Bardeen’s Fq. (31) as 


-jw 
ng® = sen(z) see J F(s) K(s) ds, ka <1 
kp 0 (3.3-4) 


where 


-jkA, ~jkA, 
K(s) =e sinh(kB, ) sinh(ks) - e sinh(kB,) sinh(ka) 


sinh(ka) 


sen(z) 


-l1,2<¢0 
€  * permittivity of free space (farad/m) 


/? 24/2 


ie - sgen(z) 2 ira cos8 - a | 


A, - 48, = (p° + [-Ja +2 sgn(z)]°}° 


Ay - 1B, = (p? + [-jis +z aanCes poh = tr” - sgen(z) 2irs cos@ - gle 


s = dummy variable with the dimension of length 
F(s) = excitation function related to the radial electric 


near-field intensity ge which is tncident on the 


groundplane. 
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The second term of Fq. (3.3-1) is therefore given by 


-Jw * 
a = sen(z) — lim F(s)K(s) ds, ka 


IA 


v0 


-jw - 

= spen(z) 2 f [ 14 r(s)| K(s) ds, ka ¢ 1 
p v0 ie 

0 (3.3-5) 


The excitation function F(s) is given by Bardeen’s Eq. (35) as 
fo. @) 
2 2.1/2 
F(s) Se se 2 Pp h(p) sinf[s (p k ) ] dp (3.3-6) 
wt 2 V2 1/2 
0 (P ) 
where 


p = dummy variable with the dimensions of (length)~! 


h(p) = function of the radial electric near-field intensity 
pfe)) 
p 


2n0 which is incident on the groundplane. 


The function h(p) is given by Bardeen’s Fq. (33) as 


ore) 
a Po lax -f h(p) J, (pp) = dp (3.3-7) 
0 


where J, (x) = Bessel function of the first kind. 


The electric field itntensity au eer incident on the groundplane 
is given by 


Bl = Ge Lee x A) (3.348) 


where 


aCe) = magnetic vector potential generated by the monopole 
source points 
H, = permeability of free space (henry/m) 
@_ = unit vector in the cylindrical radial direction 


p 


The magnetic vector potential ae) for the element base at an 


arbitrary height v is given Bo 
h 
ate) 2 cari f 1(8) (1/7Q) exp(-JkPQ) dé (3.349) 
0 
where 
/~ 
1 = unit vector along the z axis 


PQ = distance from source point Q to field point P 
mp? + (2 - & - v)y1/2 , acxt, v7, €) OL, 9, 6). 
1(€) = monopole current distribution = [1(0)/sin(kh)] sinfk(h-é)]. 


In order to facilitate the evaluation of h(p) in Fq. (3.3-7), it 
is convenient to express the factor (1/P?Q) exp(-4kP0) in terms of 
Bessel functions. Using Sommerfeld’s formula(2!) and the dummy 


variable p introduced in Fq. (3.3-6), 


[o.@) 
en 1kPQ J,(p ) exp | -lz2-é- vl vp’- «| p dp, 
* a 5 2y 172 
—— eo) 
- (7/2) < arg kad < (7/2), -7K< arg k < 0 
(3.3-19) 


Substituting Eq. (3.3-10) into Eq. (3.3-9), 


h foe) 
=e) M, 1(0) J (pp) exp |- Iz -é- v Vp? = 7 
a°'s Garataceny 2 f48 etatkCé — 6) 231? pops 
0 0 Cpe ek) 
- (7/2) < arg ia - me < (7/2), -< arg k < 0 (3.3-11) 


(e) 
p 


z<é+v. Substituting Fq. (3.3-11) for such a case into Eq. (3.3-8) 
and noting that 


Since we are interested in evaluating EF 220° consider the case 


~ x oJ. (pp) 
ae) “= A and = = pJ,(pp), 
2 
ne) : 1 rs 
p z=0 “Jenne, dzdP z=0 
CO 
7 1(0) p's, (pp) exe 2-v) Vp" - «| R(p) dp 
= 547 WE sin(kh) Pp in 
0 Zz 
where (3.3-12) 
h 
B(P) -f sin{k(h - &)] exp | -§ on «? | ag 
0 
_ _k exp { -hvp” =k) + vp" ~ k? sin(kh) - k costikh) 
2 
P (3.3-13) 


Substituting Eq. (3.3-13) into Fq. (3.3-12) and letting z=, 


fora) 

ze 2 “2 = 2/3 we 
p(e = (0) dp J, (pp) e e 

P |o20 ~j4Twe, sin(kh) 


0 


+ vy = Ke sin(kh) - k cos(kh) | (3.3-14) 


a 
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Comparison of Eq. (3.3-14) with Eq. (3.3-7) yields 


| har (h + v) Vp" - «? 


1(0) 
janwe. sin(kh) 


pe cae ee 
+ Vp" - Ke? sin(kh) e ae k - k cos(kh) e ce : 


(3.3-15) 


h(p) = 


Substituting Eq. (3.3-15) into Eq. (3.3-6) 


F(s) = a [hae aint) I, - cos(kh) 14] 
(3.3-16) 
where 
ee 1(0) 


j AWE sin(kh) 


i [2 2 
lL = p e fh ey NDE ae sin(s a - K2) dp 
1 ( 2 ) 1/2 


i) Peek 
[o¢) 
2 2 
1, = p sin(s me 2) e” Se erie dp 
(s) 
co /? 
2 2 -v™” - «2 
I. = p sin(s Yp -k ) e dp 
3 


f (>°- x’) 1/2 


In integrals Ty, 19, I,, we introduce a change in variable p to y given hy 


1/2 


y= (p - k’) - (1/2) < arg y & (7/2) (3.3-17) 


where the condition on y follows from the conditions of Eq. (3.3-10). 
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From Integral 2.663.1 of Ref. [22], it follows that I, ts given by 


io.) 


Il. = eon + v)y eid(ay)-dy = eat 7 yy [-(h + v) sin(sy) 
1 ae a ee PORE | 
jk (h + v)"+ 8 
= -k(h + v) 
~ s cos(sy)]] = —t—-; —5- {-3(h + v) sinh(ks) - s cosh(ks)] 


(v + h)” + 8 


jk (3.3-18) 


From Integral 2.667.5 of Ref. [22], it follows that I, is given by 


00 
_ SONY, ? 2 
I, = y sin(sy) e be dy = -vy + ea = sin(sy) 
2 2 2 2 2 
vts vts 


jk 
9 ore) 
- € + ur) coscer| 
stv tk 
_ ~jkv 2). 2 7 
_ e 5 i sinh(ks) - jks cosh(ks) - j(v s ) sinh(ks) - 2 vs cosh(ks) 
stv stv 

(3.3-19) 

From Integral 2.663.1 of Ref. [22] it follows that I, is given by 

fore) 
I, = f sin(sy) e dy = wor] wats) senate) | 
jk stv tk 
(3.3-20) 
-jkv - j v sinh(ks) - s cosh(ks) 
=z -e ——-:_ —___—- 
2, 2 
stv 
Consider now the limits of I> Ty, and IT, 48 v0. 
e ikh 
lim I, = —>z 5 _ [-3 h sinh(ks) - 8 cosh(ks)] “3 
seer 1 h2+ 32 (3.3-21) 
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lim I, = -7kd(s) sinh(ks) + ik coshtka) 


v—0 : 
2,2 
+ jv 6° (s) sinh(ks) - —isinhtks) -7§6°(s) cosh(ks) 
3 
(3.3-22) 
lim I, = 476(s) sinh(ks) + —coshtks) (3.3-23) 
v0 
ehapecre? 
é(x) = Dirac delta function = lim E az | 
re 2 2 
v0 v+x 
b 1, a< Xo <b 
[ b0-«,) dx = 1/2, as x <bora<x=b 
a 0, x €aorx >b 
fo) fo) 
Rye (800) .g. tae 3 [ 1 v: |- lim - 1 2 vx 
dx v0 Ox T yet x2 v—0 v (v7+ x”) 2 


In Eqs. (3.3-22) and (3.3-23), those terms containing the product 
6(s) sinh(ks) may be set equal to zero since from Fq. (8. 2-12) of Ref. 
{23], x 6(x) = 0. 


The lim F(s) is found by substituting Fqs. (3.3-21)-(3.3-23) into Fq. 


vo 


(3.3-16). 
Accordingly, 
-4jkh 
2Ck e [-i h sinh(ks) - s cosh(ks)] 
pe ga 
v—0 h+s 
. 2C sin(kh) -j k cosh (ks) _ 4 sinh (ks) 478 (8) cosh (ks 
1 s 2 
s 
+ —2C Kcoe(kh) costi(ks) (3.3-24) 
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The magnetic field intensity generated by the groundplane current, for 
the case v0, is found by substituting Eq. (3.3-24) into Eq. (3.3-5). The 
total magnetic field is explicitly determined by numerical evaluation of Fqs. 
(3.3-1) - (3.3-5). 


In the far-field, the factor K(s) and the parameters Ay, Ay, By, Ro, in 


Eqs. (3.3-4) and (3.3-5) reduce to a simpler form. 


In the far-field when r—aw, then a << r and s << r. For these 


conditions, 


A, - 5B, Xr" - sgn(z) 2jra cos 8 - a?cos2a}!/2 = r- sgn(z) ja cosé@, 
ar (3.3-25) 
Ay - jBy RE - sgn(z) 2jra cos 6 - s?cos-a}!/2 = r- sgn(z) is cosé@, 


s <r (3.3-26) 


Equating real and imaginary parts in Eqs. (3.3-25) and (3.3-26) 


A, = Ay = 03 ar, sr (3.3-27) 
B, = sgn(z) a cos 8, a <<p (3.3-28) 
Bo = sgn(z) s cos 8, s <<p (3.3-29) 


Substituting Eqs. (3.3-27), (3.3-28), (3.3-29) into Fqs. (3.3-4) and (3.3-5), 
the far-field magnetic field intensity due to the groundpJane current, for 


the case v—0, is given by 


~jwe exp(- jkr) 2 


lim wi = Se Of tie r¢s)] g(s) ds; far-field, ka ¢ 1 
q: hve (3,3-30) 
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where 


sinh(ka cos8) sinh(ks) - sinh(ks cos@) sinh(ka) 


g(s) = sinh(ka) 


and lim F(s) is given by Fq. (3.3-24) 
v0 


Eq. (3.3-30) can be reduced further by utilizing the properties of the 
Dirac delta function given by Fq. (B.2-7) of Ref. [23] and the relation 


2 1 
f 6°(x) f(x) dx = - “= £7(0' 


0 
It follows that Eq. (3.3-30) reduces to 
(2) 2A | a e Jkh 
lim H.® = —— ds g(s){—~——s—__ [-41 h sinh(ks) - s cosh(ks)} 
) T 2 2 
v—0 h +s 
0 
+ sin(kh) -4 cosh(ks) ~ 4 #ion(kay 4 cos(kh) cosh(ks) 
s 2 s 
ks 
sinh(ka cos @) - cos 4 sinh(ka) 
- OTRAS, COS COS SER Re? - 1 
a,stnceh)| sinh(ka) ; far-field, ka < 


(3.3-31) 


where A and g(s) are defined in Eq. (3.3-3) and Fq. (3.3-30), 
respectively. 


The total far-field magnetic field intensity generated by the 
element and groundplane currents, for the case v-0, is found by 
substituting Eq. (3.3-3) and Fq. (3.3.-31) into Fq. (3.3-1). The 
(total) { 


resulting expression for lim H 
v0 


s in a form suitable for 


computer evaluation. 
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The radiated power density s(8), directive gain d(8), and 
radiation resistance R are found from Eqs. (2.2-4) and (2.2-7). These 
quantities are computed in program “BARDEEN” in Appendix Bl. Computer 
printouts of the relative power radiation pattern s(8)/s, 4, radiation 
resistance R, and the directive gain d(®) are given in tables (Al-1 - 
Al-13), (Al-14 - A1-25), (A1-26 - A1-37), (A1-38 - Al-49) for 
normalized element lengths h/A = 0.25, 0.1, 0.040, 0.01, respectively, 


and normalized groundplane radii = 0, 0.25, 0.50, ... 3.0 wavenumbers. 


The radiation patterns have no appreciable change in shape for 
groundplane radii 0 < kaX<1.75 and resemble that of a dipole in free 
space with peak gains approximately in the direction of the horizon 
and with directivities less than that of a dipole whose total element 
length is twice that of the monopole element length. For a 
quarterwave monopole element, the directive gain on the horizon 


deceases from 1.88 dBi for ka = 0 to 1.23 dBi for ka = 1.75. 


In table 5 of Section 3.9, the radiation reststances obtained by 
the integral equation method are compared with those obtained by the 
method of moments for a quarterwave element and groundplane radii 
0 < ka ¢ 3.0. The values, obtained by the integral equation method, 
differ from those obtained by the method of moments’ by less than 
1% for 0 < ka ¢ 1.75 and by less than 10% for 0 < ka ¢ 2.75. These 
results suggest that the integral equation method is accurate for 
groundplane radii 0 < ka < 1.75 and is useful for 0 < ka ¢ 2.75. It 
is suspected that the best available results are obtained by the 
integral equation method for 0 < ka ¢ 1.0. 
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3.4 Method of Moments, 0 < ka <14 


In Section 3.3 it was shown that the integral 
equation method, for determining the fields generated by the 
groundplane current when the current distribution on the monopole 
element is initially specified, is accurate only for groundplane radii 
less than approximately 1.75 wavenumbers. In Section 
3.5 it will be shown that the algorithms utilized in the oblate 
spheroidal wave function methn4 are accurate only for groundplane 
radii no smaller than 3.0 wavenumbers. The question arises: What 
method is accurate over a range of groundplane radii which includes 
the region 1.75 < ka < 3.0? 


We have found that Richmond’s method of moments (2) is the only 
method in the present literature to be accurate over a range of 
groundplane radii which includes the region 1.75 < ka < 3.0. Although 
this method is primarily intended for use when the current 
distribution on the monopole element is initially unknown, this method 
is also applicable when the element is specified to have a sinusoidal 


current distribution. 


This method is discussed in Section 4.2. In this method the 
element is subdivided into N equal segments and the groundplane is 
subdivided into M concentric annular zones of equal width. The 
unknown current distributions on the element and groundplane are 
expanded as a series of N + M overlapping sinusoidal dipole modes 
(sinusoidal-Galerkin method) each with an unknown current 
distribution. The N + M currents are determined by inversion of a 


(N + M) x (N + M) matrix. The numbers of subdivisions, N and M, are 


limited by the cost of computation time and by the precision of the 
computer. The accuracy of the solution can be decreased appreciably 
if either N or M is too small or too large. The method of moments 
converges to a solution when an increase or decrease of unity in the 
value of N or M does not appreciably alter the solution for input 


impedance. 


A particularly useful property of the sinusoidal - Galerkin method 
is the sinusoidal current distribution which is imposed on the element 


by setting N=l1. 


The input impedance and radiation pattern of thin quarterwave 
elements (b/A= 10-6, h/A = 0.25), for groundplane radii 0 < ka < 14 
wave numbers, were determined by Richmond’s method of moments by 
utilizing MITRE programs “RICHMD1” and "RICHMD2” which are discussed 
in Section 4.2 and listed in Appendix B5. With N=l, convergent 


solutions were obtained for values of M given by: 


M = 3, ka = 0.25 

M=7, ka = 0.50 

M = 16, ka = 0.75, 1.0 ..... 8.5 
M=3 ka, ka = 8.75, 9.0, ..... 14.0. 


Computer printouts of the input admittance, groundplane current 
distribution, radiation resistance, directive gain, and radiation 


patterns are given in Appendix AS. 
The input resistance and reactance -- as determined by the method 


of moments or methods which give similar results -- are plotted for 


thin quarterwave elements on groundplanes of radii 0 < ka < 14 in 
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figures 9 and 10 of Section 3.9. The directive gain on the horizon, 
peak directivity, and elevation angle of the peak directivity given by 
the computer printouts of Appendix A5 are plotted in figures 11 - 13 
of Section 3.9. 


The input impedance of thin elements with a sinusoidal current 
distribution were also determined for element lengths h/A = 0.1 and 
0.025 and groundplane radii 9 < ka < 8.0 wavenumbers by utilizing 
program RICHMD1 with N=l. The radiation resistance of these elements 
is compared in figure 14 of Section 3.9 with those of a quarterwave 
element. In figure 14 the radiation resistance is normalized to the 


value of radiation resistance of each element for ka = 0. 


Values of the radiation resistance for quarterwave elements are 
determined in Appendix A5 by matrix inversion (program RICHMD1) and 
also by the far-field radiation pattern (program RICHMD2). The values 
determined by both methods differ by less than 1% for small 
groundplane radii and differ by less than 3% for the larger 


groundplane radii. 


The values of radiation resistance obtained from program RICHMD2 
are compared in table 5 of Section 3.9 with the values obtained by the 
{integral equation method and the oblate spheroidal wave function 
method. Richmond’s method of moments gives useful results over the 
entire range 0 < ka < 14 and gives good agreement with the integral 
equation method for 0 < ka ¢ 1.75 and with the oblate spheroidal wave 
function method for 3.0 < ka < 6.5. Whereas the method of moments 
gives useful results in the regions 1.75 < ka < 3.0 and 6.5 < ka < 14, 
the other methods fail In these regions. For ka > 14, Richmond’s 
method of moments Is not as useful because of increased computation 


time and decreased accuracy. It is suspected that for thin monopole 
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elements the best available results are obtained by the method of 
moments for 1.25 < ka < 2.75 and 6.75 < ka ¢ 14. For relatively thick 
monopole elements, Richmond’s method of moments gives the best 
available results for groundplane radii 0 < ka <¢ 14 as discussed in 


Section 4.2. 
3.5 Oblate Spheroidal Wave Functions, 3.0 < ka < 6.5 


Oblate spheroidal coordinates (see table 1 in Section 2.1) are 
particularly convenient for handling the boundary conditions of the 
magnetic field intensity on the groundplane. The requirement for 
constant tangential magnetic field intensity across the upper and 
lower hemispheres at the groundplane interface may be specified at all 
points of the groundplane disc by a boundary condition at only the 


oblate “radial” coordinate & = 0. 


For groundplane radii of the same order of magnitude as the 
excitation wavelength, Leitner and spence’?)~ (5) utilized oblate 
spheroidal wave functions to determine the groundplane current induced 
by a thin quarterwave element with a sinusoidal current distribution. 
Leitner and Spence give numerical values of the groundplane current 
distribution, radiation resistance, and far-field power density (at 
constant element base current) for groundplane radii ka = 3, 4, 5, (42 
wavenumbers. The complex current distributions on both the top and 


bottom faces of the groundplane are reported. 


In this section, we report the results of a computer program 
“MONOPL," based on the theory of Leitner and Spence, which calculates 
the directivity pattern and radiation resistance for groundplane radii 
including the cases ka = 3, 4, 5, (42 wavenumbers. Our results are 
consistent with (but not identical to) the results reported by Leitner 
and Spence. The form of the solution, corrections, regions of 
calculation validity, accuracy, and numerical results of program 


MONOPL are discussed in the remainder of this section. 
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Form of Solution 


For a quarterwave monopole (h =A/4), the far-field power density 
s°(7) and radiation resistance R, expressed in oblate spheroidal wave 


functions, are given with reference to figure 1 of Section 2.1, by) 


foe) 
s(n) = [12(0)/8r?r?) uife.)/? e4 | SB, wy, CDI? (3.51) 
f=1 


‘et ae eee 2 
R= (1/27) (u,/€,) /2 ¢ Ny, By!’ obms (3.5-2) 
f=1 


where 
= wave impedance in free-space (ohms) 
€ = ka = 27a/X = groundplane radius (wavenumbers) 


uy, (7) = angular oblate spheroidal wave function of 


order one and degree f 


Ney -f ug {DAN = norm of Uy, (7) 
= 


oY) = + g2yl/2 sete) 


(3) 


= f21)271 2 
By = (“1 gh yy CT/ZED$ONp 9p) Mg 09) 


The reader ts referred to reference [3] for a definition of the oblate 


(3). 3 
spheroidal wave functions @),, q)1, vp, (é), and 9%16$)> 
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Substitution of Eqs. (3.5-1) and (3.5-2) into Eq. (2.2-4) yields the 
directivity d°(7) given by 


[o) io) 
am) = 215 Bul? / Nyy Iw, I? (3.5-3) 
fel £=1 


Corrections 


Problems were encountered in generating the eigenvalues of the 
oblate spheroidal wave functions. These problems were narrowed down 
to the continued fractions that generate the elgenvalues. These 
continued fractions were derived from the recurrence relations which 
in turn come from the differential equation. A discrepancy was 
discovered between our derived continued fractions and the published 
continued fractions [3] . The corrected continued fractions appear 


below. 
g-m even 


-W4m) (a-m) €2 


em” 2(2e=1)4 SOS D)_(e-m~2) «2 ; 
y 4(2£-3)+¥. +  (L4+m-4) (g-m-4) € 
£m 6(2i-5)+y + 
£m 
: : 
vs (é+m+2) (g-m+2) € 2 
2(2£+3)-y, + _Cétmt4) Geom 4) €O 9 
Ee B(W+5)- y + _@mr6) (2-m+6) —€* 
ay 6(2i+7)-y. OF 
£m 
(3.5- 
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£- m odd 


y 2 ribs) (em)? ; 
£m 2(2k-1)+ y + (£+m-3) (£-m-3) € > 
fh 4(28-3)+ 7 + (g+m—5) (g-m-5) €° 
6(20-5)+ ¥, — + 

(g+mtl) (£-mt1) €2 ; 
2(24+3)-Y + (£+m+3) (£-mt+3) € : 9 
k 4(2E+5)- 9, + ES) (u-mt5)_ 
6(264+7)-y + 

(3.5-5) 


The eigenvalues can also be expressed in terms of a series expansion, 
which has the form 


lo. 4) 

gm 2k = 

Ym 7 > fe OE (3.5-6) 
k=1 


The first two terms in this expansion were checked against the continued 
fractions and were found to agree. Thts is important, since the continued 
fraction method by which values of the eigenvalues are obtained depends upon 


the accuracy of the roots 7, in equation (3.5-6). 


Regions of Calculation Validity 


A lower bound on the value of € = 27a/, for which the calculations are 


valid depends upon the following equations [see Fqs. (34) and (48) of 
Ref. 3]: 


[o.@) 
Tie) = (14421 /2 2 ' as en for f= 7/2€ (3.5-7) 


n=c, 1 
and 


21/2 Al 
B,(E) = (14é*) ane a 


n=0,1 


5 for t= 1/2€ (3.5-8) 


where the prime indicates summation over alternate n, starting with n=0 if 
(£- 1) is even, with n=l if (£- 1) is odd. 


Eq. (3.5-7) is the expression for the first order radial spheroidal 
wave function of the first kind. Eq. (3.5-8) is used in another 
expression to obtain the first order radial spheroidal wave function 
of the second kind. As € becomes small, § becomes large, so more 
terms are needed in Eqs. (3.5-7) and (3.5-8). In theory, Eq. (3.5-7) 
and (3.5-8) converge for all real values of = 7/2€. Computationally, 
however, because of the finite accuracy of the computer, Eqs. (3.5-7) 
and (3.5-8) will not converge for all real values of 

§=7/2€. To be on the safe side, we can restrict & so that é is less 
than unity. With this assumption, we get the computational constraint 
€=7/2€ < 1 which implies a lower bound on € given by € >7/2%1.57. 


An even tighter lower bound on € is obtained by observing what 
happens in the algorithm that is used to obtain the eigenvalues. This 
algorithm is not well behaved for values of € less than € = 2.5. 
Therefore € is lower bounded by € > 2.5. For accurate values of 
radiation resistance, € is lower bounded by € > 3.0 (see table B of 
Section 3.9). 
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An upper bound on € is obtained by observing what happens in the 
continued fraction algorithm that is used to obtain the eigenvalues. 
For values of € greater than 6.5, the series expansion for the 
eigenvalues, given by Eq. (3.5-6), does not give an accurate enough 
answer using only the first four terms. The resulting values for the 
eigenvalue are far enough from the correct eigenvalues to cause the 
continued fraction algorithm of Eqs. (3.5-4) and (3.5-5) to converge 
to a root which is not the eigenvalue. As a consequence, erroneous 
values of radiation resistance can be obtained for €> 6.5 (see table 5 


of Section 3.9). 


Consequently, the range of € for which Leitner and Spence’s method 
of oblate spheroidal wave functions is useful is 2.5 < € < 6.5. It is 
suspected that the best available results are obtained by the method 
of oblate spheroidal wave functions for 3.0 < € < 6.5. 


Accuracy 


There are many equations involved in the calculation of 
directivity and radiation resistance. Some of these equations involve 
series expansions. When we varied the number of terms in these 
series, the radiation resistance was found to vary only in the fifth 


or sixth significant figure. 


Another problem that was mentioned previously was the accuracy by 
which Eq. (3.5-6) computes the eigenvalues. Because of computational 
reasons, only the first four terms in Eq. (3.5-6) were used. Because 
of this, a raw eigenvalue is computed using Eq. (3.5-6), which is then 
used as an initial guess in the continued fractions. The continued 
fractions have many roots in ae The number of roots is dependent 
upon the number of fraction terms used in the continued fraction. 

Only one of these roots however can be the eigenvalue. if the raw 
eigenvalue is far enough fron the correct eigenvalue then the 
continued fraction will converge to a root which is not an eigenvalue. 


This will produce wrong results. 
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The eigenvalues that were computed were checked against published 
values**”? . The computed values were found to be within the 


percentage error of the published values. 


The computed values of directivity on the horizon and radiation 
resistance, for different values of groundplane extent, were found to 
agree with the relationship given by Eq. (2.3-3) to at least five 
places after the decimal point for the free space wave impedance 

1 = 376.73037 ohms. 


Numerical Results 


The far-field power density s°(”7), radiation resistance R, and 
directive gain d°(7) of quarterwave elements given by Eqs. (3.5-1), 
(3.5-2) and (3.5-3), respectively, were numerically evaluated by 
program MONOPL written in FORTRAN 77 language for use on a DEC 
PDP-11/70 computer. The program listing is given in Appendix B3. 
Numerical values were obtained for the cases 2.5 < ka < (42. 


The computed eigenvalues are given in table 2. Computer printouts 
of the directivity patterns are given in Appendix A3 for ka = 3, 4, 
5, /42. The patterns are plotted in figure 8 of Section 3.9 as polar 
graphs on the same linear scale. In these plots the total radiated 


power is held constant. 


The directional gain on the horizon, radiation resistance and peak 


directivity are summarized in table 3. 
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_— 


1 
2 
3 
4 
5 
6 
7 
8 
9 
0 


Se ee ae ee ae ee a OL oD 


Table 2. 


- 7305791912 
- 8215414322 
0070741733 
~ 7772522576 
- 6896340053 
-6368101196 
- 6032991570 
- 5807102835 
-5647761329 
-5531234272 


Eigenvalues 


1] 


WwemOmDmDDOaWWO Oo 


71 For 


- 0907992567 
- 1301387077 
- 4059031099 
- 6287184908 
- 4462229270 
- 3210439564 
- 2427831207 
- 1899480807 
- 1526212369 
- 1252813999 


2n a/A = 3,4,5, 7 42 


056611686 
-872162474 
- 750367134 
-694959245 
- 439315734 
- 161309568 
-999910517 
-890952874 
-814155253 
- 757940027 


- 687439619 
- 7382 30003 
127723151 
- 237303490 
- 585994650 
- 624607030 
- 234563635 
- 960565094 
- 770549914 
-632114595 


€9 


Table 3 


Directive Gain on Horizon and Radiation Resistance of 
Quarterwave Element, 2.5 < 2ma/A< /42. 


Peak Directivity 


Leitner* (dBi) 


Directive Gain 
on Horizon 


Radiation 
Resistance, R 


Groundplane 
radius, 


€= 2ma/Xd 


d(7/2) D(7/2) = 


(Wavenumbers ) 


(numeric) 


10 log, )4(77/2) 
(dBi) 


Program 
MONOPL 


re 1.09745 0.40383 27.31731 
3.0 0.79995 -0 96938 37 47646 2.52250 
3.4 0.66233 -1.78923 4526299 
3.5 0.65086 -1 86515 46 .06128 
3.6 0.64812 -1.88346 46.25589 
3.7 0.65316 -1.84980 4589866 
4.0 0.70254 -1.53326 42.67237 3.89943 
4.5 0.82743 -0.82271 3623191 
5.0 0.92059 -0 35934 32.56526 3.37175 
5.25 0.93879 0.27430 31.93379 
5.5 0.93294 -0.30146 32.13410 
6.0 0.85088 -0.70134 35 .23337 
J42 0.74637 1.27045 40.16661 2.55213 


* Research Report No. EM-19, New York University, Washington Square College, April 1950. 


3.6 Scalar Diffraction Theory, Geometric Theory of Diffraction, 
6.5 < ka < 


For an element on a groundplane of sufficiently large radius, 
Tang‘) utilized a scalar theory of diffraction to calculate the 
far-field elevation pattern. For elevation angles near the horizon 
(97/2) the far-field magnetic field intensity is determined by 
linear extrapolation to the result for the element itself given by 
Eq. (3.2-10) with 9= 7/2 rad. Tang“s method for the radiation pattern 
is more accurate than that obtained by the variational method of 
Section 3.7 because it includes an additional term in the expansion 
for the total magnetic field intensity. Since the variational method 
is useful for groundplane radii as small as ka=30 wavenumbers, Tang’s 
method should be useful for even smaller groundplane radii provided 
ka >> 1. 


The geometric theory of diffraction (GTD) is another method which 
is applicable for sufficiently large groundplane radii. in GTD, the 
fields are determined by ray optics (an incident ray plus a reflected 
ray) and diffraction by the edge. However, the effect of edge 
diffraction {s only approximated because in GTD the edge is treated 
point by point as though it were a straight knife edge of infinite 
extent. For this reason, GTD may be applied to an element at the 
center of a circular groundplane only when the groundplane is of 
sufficiently large radius. The method of GTD is reviewed by 


(27) 


Balanis who also gives a computer program for calculating the 


diffraction coefficient. 


In Section 4.3, the method of moments combined with GTD yields 
results for input impedance which are useful for groundplanes of radil 
ka 2 6 wavenumbers and are accurate for ka > 8. Therefore, when the 
element current distribution ls constrained to be sinusoidal, the 
method of GTD ls expected to give useful results for groundplane radii 
ka > 6.5 wavenumbers over which range the nethod of oblate spheroidal 


wave functions does not give useful results. 
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3.7 Variational Method, 30 < ka < 


For an element on a groundplane of very large radius 
(ka > 30 wavenumbers), differences in the input impedance and 
radiation pattern from that for a groundplane of infinite extent may 


be determined by utilizing a varfational method of Storer(7) »(8) 
With reference to figure 1 of Section 2.1 and for a sinusoidal 


current distribution on the element given by Fq. (3.1-1), the input 
impedance difference is given by Fq. (20) of [7] as 


1 + 


(uma)! 
ka > 30 (3.7-1) 


| a ae -j" exp( j2ka) | 1- eeege 
in oo aka sin(kh) 


exp[i(2ka + 3971/4 1 


where 


Zan = input impedance for an element on a groundplane of 
radius ka wavenumbers (ohms) 

Z_.= input impedance for the same element on a groundplane 
of radius ka =OOwavenumbers (ohms) 


1} = wave impedance in free space = 376.73037 ohms 


Since (anka)!/? >> 1 for ka > 30, the input resistance difference 


Rin 7 Roo and input reactance difference X,. - Xo are given 


in 
approximately by 


2 
1 - cos(kh) | sin(2ka) 

R.- R =| Senn ine) ka > 30 

in 00 sin(kh) 4Ttka (3.7-2) 
1 eesthhy. cos(?ka) 

x -x wx|—-tccostkh) pei A AS PI ka > 30 

in “oo sin(kh) | Uttka , aay tee 
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The maxima, nulls, and minima of RR, - R_ occur for values of ka given 


in 
approximately by 


(2N + > ) (7/2) (maxima) 
ka = N (71/2) (nulls) » N = 0,1,2,.... 
(2N + = )(71/2) (minima) 


(3.7-4) 


The input impedance Z,_ -—- Z, given by Eq. (3.7-1) ls calculated in 


in 
computer program “MONOSTOR” whose listing is given in Appendix B-4. 


For very thin quarterwave elements, Ry* 36.54 ohms and 
: re 21.26 ohms (see Section 3.8). The numeric directivity on the 


in by 


Eq. (2.3-2). Computer printouts of the maxima and minima of Rin = Ry 
X 


horizon, d(7/2), is related to the radiation resistance R = R 


ine Xp are given in Appendix A4. For very thin quarterwave 


elements, Appendix A5 also gives computer printouts of Rea x 
d(7/2), and p(7t/2) = 10 log 10 d(7/2). 


in’ 


Differences in the far-field radiation pattern from that for a 
groundplane of infinite extent are given in [8]. For the waveform of 
Eq. (2.1-2) and an element sinusoidal current distribution I(z) (see 
figure 1 of Section 2.1), the difference in the far-field magnetic 
field intensity is given by Eq. (6) of [8] as 
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4 .1(0) exp(-jkr) [1-cos(kh)] | 


sin(kh) 27r °- 27/2 
27 
r a = 
if exp[-jka(1-sin® cose) cosh a] , ka > 30 
0 (l-sin® cosd) (3.7-5) 
where 
Ng = far-field magnetic field intensity for an element on a 
groundplane of radius ka wavenumbers (amp/m) 
Ho! paw co = far-field magnetic field intensity for an element on a 


groundplane of infinite extent (amp/m) 
#1, 0<0<7/2 
sgn 8 = 0 


-1,7/2<e<7 


The magnetic field intensity Hg! ax cots given by (see Section 3.8) 


1(0) exp(-jkr) cos(kh cos8) - cos(kh) } 0<e<7/? 
sin(kh) 27r sin @ 


0,7/2<0<7 


Hg! kaso 


(3.7-6) 


The far-field electric field intensity Eg = 1H where Hg is given by 


Eq (3.7-5) and ") is the free-space wave impedance. 


One of the distinguishing features of the far-field radiation 
pattern for groundplane of large radii is the occurrence of a 


fine-structured lobing pattern, Ng 7 Hy| superimposed on the 


=0? 
pattern for a groundplane of infinite nae The lohing pattern, 
given by Eq. (3.7-5), is symmetrical about the horizon (expect for a 
phase shift of 7 rad) with the most prominent lobes near the zenith and 
nadir directions. The a" maximum of the lobing pattern in each 
quardrant decreases with increasing values of n where n=l corresponds 


to the lobe nearest the zenith or nadir direction. 
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The lobing structure is characterized in [8] by the elevation 


th maxima, the angular separation A@ between Jobes, 


angle 8, of the n 
the angle ae within which all prominent maxima occur, and the number 
N of prominent lobes. The elevation angle 84 in radians is given 


approximately by 


0.59 (7/ka), n=1 


© 
v 


, ka > 30 (3.7-7) 
(n + 0.75)(7/ka), n = 2, 3, .2..N 


The angular separation A9 in radians is given approximately by 
1.16 (7/ka), n=1 
Ae = » ka > 30 (3.7-8) 


Wika, n = 2, 3, «+N 


Prominent lobes are defined in [8] as lobes whose maxima are less than 


one-fourth the amplitude for that of a groundplane of infinite extent. 


The angle De in radians is given approximately by 


1 - cos(kh) ie 


1/3 
8 = 1.87 (7/ka) [ eR ka > 30 
max kh sin(kh) ~ (3.7-9) 
The number N of prominent lobes is determined by 
8 = 6, + (A®)__, + (N-2)(A@) , ka> 30 
max 1 n=1 n*1 (3.7-10) 


where @,,48, and Buax ate given by Faqs. (3.7-7) - (3.7-9), respectively. 
Solving for N, 
-2/3 
N= [0 .,/(7/ka)] - 0.25 = 1.87(77/ka) -0.25, ka > 30 


(3.7-11) 
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3.8 Method of Images, ka = © 


For the idealized case of a monopole element mounted on a groundplane 
of infinite extent and of infinite conductivity, the monopole antenna may 
be modelled by the method of images as a dipole with one-half the input 


impedance and double the peak numeric directivity of the dipote 6°”? 


The 
infinite groundplane prevents monopole radiation into the hemisphere helow 
the groundplane but allows a radiation pattern identical to that of the 


dipole in the upper hemisphere. 


In this section, it is assumed that the current has a waveform given by 
Eq. (2.1-2) and a current distribution T(z) on the element and its image 
(see figure 7) given by 


[1(0)/sin(kh)] sin{k(h-z)], 0 < z < h (element) 
I(z) = (3.R8-1) 
[1(0)/sin(kh)] sin{k(ht+z)], - h <2z <0 (image) 


The near-fields, far-fields, and input impedance of an electrically-thin 
element on a groundplane of infinite-extent are summarized in the remainder 


of this section. 


Near-fields 


The exact magnetic field intensity H = Ugh and electric field 
intensity E = a E, + aE, for an infinitely thin element are given with 


reference to Fqs. (7-15) and (7-17) of [25], as 


4Tp sin(kh) 


4 1(0) aa + exp(—4KR, ) - 2 cos(kh) cxo(—e)] 


z 2 0 


0,z<o0 
(3.8-2) 
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Figure 7. 


GROUNDPLANE OF 
INFINITE EXTENT 


P(r,0, d) 


Method of Images for a Monopole Element at the Center 
of a Circular Groundplane of Infinite Extent 
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(z+h) exp(-4kR, ) 


into) | (z-h) exp(—1kR, ) . 
“GI p sin(kh) R, 
- 22 cos(kh) exp(-ikr) a, (3.8-3) 
r i 
0,z<0 
_-3n 10) _ exp(—4kR, ) P exp(— 1kR,, ) 
Gir sin(kh) oe 
- 2 cos(kh) exp(-4kr) 50 C384) 
r BT eee 
0, z2<0 
where 
R= 10 + (en)? 
2.1/2 


Ut 


R, = [D> + (2th)”] 


The magnetic field intensity at the top and bottom surfaces of the 
groundplane are given by Eq. (3.8-2) as 


eg AOD, és 2 = 
Hg! =o, rapes, [enne 4kR, ) cos(kh) exp( Hp) 
top of groundplane (3.8-5) 
Hy! = = 0, bottom of groundplane (3.8-6) 
where 
R, ~ (p? + n2y1/? 


The incremental rsroundplane current AI(P) contained within a differential 
azimuthal angle db may be de. ermined from Fqs. (3.8-5), (3.8-6), and 
Ampere’s circuital law applied to a closed path along the top and hottom 
surfaces of the groundplane along arc lengths within dh. Since the 
groundplane current I(P) is defined to be positive in the positive Pdirection 


(see figure 1 of Section 2.1), the path is taken in the clockwise direction. 
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Accordingly, 


AICP) = (ug re si “6 “Mnottom? P ve [-Mg! ene y Hg! pen_1P a6 


(3.8-7) 
The total groundplane current T(P) is given by 
27 
1(p) = f AI(p) dd = 27O[-Hyl gy + Aylin! 
0 
j I(9) feos(kh) exp(-jkp) - exp(—4kR, ) 1 (3.8-8) 
sin(kh) 


Eq. (3.8-8) agrees with that given by Eq. (33) of [7]. 


Far-fields 


For the far-field conditions given by Fqs. (3.2-7) and (3.2-8), 
== A — 
Eqs. (3.8-2) - (3.8-4) reduce to the far-fields Heugha and Fuh gay, where, 
with reference to Fq. (4-62) in [25], No is given by 


j 109) exp(-jkr) cos(kh cos®@) - cos(kh) 
sin(kh) 27 sin@ 


Hy = (3.8-9) 
0, 7/2<0<7 


|. o<acm 


The time-averaged radiated power density s(@) =(7/2) 1g!” is given by 


nixcoy12 £(@) 
Br’? sin2(kh) 

s(®) = (3.8-10) 
Dy A/S Bee 


where 
2 


sind 


The direction of maximum radiated power is 0 = 7/? rad. 


£(6) = [ cos(kh cos®@) - cos(kh) | 
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The input impedance of a monopole element of length h on a groundplane of 
infinite extent is one-half that of a dipole of total length 2h in free space. 
Accordingly, the radiation resistance R of an infinitely thin element on a 
groundplane of infinite extent is given by [compare with Fqs. (4-79) and 
(4-79) in [25]] 


a2 n Cin(2kh) + 5 sin(2kh)[Si(4kh) - 2 Si(2kh)] 
uatsin® (kh) 


+ + cos(2kh)[2 Cin(2kh) - Cin(4kh)]}, =n 
(3.8-11) 


where Si(x) and Cin(x) are the sine integral and modified cosine integrals, 
respectively, which are defined following Fq. (3.2-32). Faq. (3.8-9) is also 
approximately valid for thin elements (kb << 1). [For example, compare 

Eqs. (3.2-33) and (3.2-35) for a thin monopole element in the absence of a 
groundplane. J 


Substituting Fqs. (3.8-10) and (3.8-11) into Fq. (2.2-8), the numerte 
directive gain d(@) is given by 


d(6) = Aft (3.8-17) 
where 
B = Cin(2kh) + — sin(2kh)[Si(4kh) - 2 Si(2kh)] 


2 


+ a cos(2kh)[2 Cin(2kh) - Cin(4kh) J 


For a quarterwave monopole element (kh=7/2), Eqs. (3.8-11) and (3.8-12) 


reduce to 
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R = i nn) = 36.5395 ohms; kh=7/2, b=0 
(3.8-13) 
2 
£(@) = £08 _[(mi2) 088} , kh=n/2 (3.8-14) 
sin 8 
4 . (71/2) ) 2 
aca) = + 08 [ON ?) coe8] | 2 C=" khen/2 
sin Cin(27r) 

(3.R8-15) 


The peak numeric directivity d(7/2) = 3.2818 corresponding to N(7/?) = 


10 10g, 4 d(7/2) = 5.1612 dBi. The directional gain is plotted in 
figure 8(f) of Section 3.9. 


For an electrically short monopole element (kh << 1) and 
7) 12077 ohms, Eqs. (3.8-11) and (3.8-12) reduce to 


R = 10 (kh)?, (kh) << 1 (3.8-16) 

4 
£(®) = ee sin?@, (kh) <<1 (3.8-17) 
d(®) = 3 sin @, (kh) << 1 (3.8-18) 


The peak numeric directivity d(7r/2) = 3.0 corresponding to D(7r/2) = 
10 log, 9 d(m/2) = 4.7712 dBi. 


Input Impedance 


The input resistance R is given by Eq. (3.8-11) which is exact 


in 
for an infinitely thin element and approximately correct for thin 


elements provided that the element current distribution is sinusoidal. 


The input reactance X for thin elements (kb << 1) is given by 


in 
{compare with Eqs. (7-33) and (7-30) in [25] 
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1 


ee n Si(2kh) + cos(2kh)[Si(2kh) - —— Si(4kh)] 
ip amt in? (kh) 
- sin(2kh)[1n(h/b) - Cin(2kh) + + Cin(4kh) + + Cin(kb2/h)] 
(3.8-19) 


For a quarterwave element, in the limit b--0, Eq. (3.8-19) gives an 


ta = 21.2576 ohms. In Eq. (3.8-19), X, 70 for 
element lengths slightly less than kh = (2n-1)(7/2), n=l, 2, 3, «.-, 


input reactance X 


which are approximately one-half the resonant lengths for a monopole 


element with no groundplane [see Eq. (3.2-38) ]. 
3.9 Summary of Results 


In Section 3 a sinusoidal current distribution is assumed on the 
monopole element. Although such a current distribution is never 
exactly realized even for an infinitely thin element (see Section 3.1), 
it is a useful approximation for sufficiently thin elements. For 
example, for a quarterwave element of radius b=107°A , the input 
impedances computed by determining the actual element current 
distribution and that computed by assuming a sinusoidal current 
distribution differ by no more than 5% for groundplane radit 
ka = 6 -— 14 (compare table 7 of Section 4.5 with tables A2-24 - A2-42 
of Appendix A). The assumption of a sinusoidal element current 
distribution allows for models which are coaputationally simpler in 
determining input impedance and radiation patterns than the models 
which follow in Section 4. The results of these simpler models are 


summarized in this section. 
The electrical properties of electrically-short and quarterwave 


monopole elements on groundplanes of zero, large, and infinite extent 


are compared in table 4 with the corresponding properties of 
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9L 


PEAK DIRECTIVITY * DIRECTIVITY ON HORIZON | INPUT IMPEDANCE | 


GROUNDPLANF. ELEMENT 
ANTENNA RADIUS LENGTH ace) D(®8_) a(@ =7/2) RADIATION REACTANCE*** REFERENCE 
STRUCTURE (WAVENUMBERS ) (WAVELENGTHS ) (numBERIC) (aBP) (NUMERIC) RESISTANCE EQUATIONS 
2ma/i h/A (OHMS) CORMS ) 


20m (h/A)? 


0 3.2-19, 3.2-20 


0 19.43 3.2-15, 3.2-16 


sor? (n/a)? 3.7-1, 3.8-18, 2.3-2 


MONOPOLE >>> 1, FINITE 


36.54 3.7-1. 3.8-15. 2.3-2 


>>> 1, FINITE 


2 


407 (h/ d) 3.8-16, 3.8-18, 3.8-19 


36.54 


3.8-13, 3.8-15, 3.8-19 


DIPOLE” g077(h/ A)? C. A. alanis [15], 


73.08 pp. 108, 112, 131, 132 


IDEALIZED ISOTROPIC RADIATOR 


* DIRECTION OF PEAK DIRFCTIVITY IS @ = 77/2 RAD EXCEPT FOR MONOPOLE ELEMENTS ON GROUNDPLANES OF FINITE BUT NON-ZERO EXTENT IN WHICH 
CASE @, < m/Z RAD. Pp 


** VERY THIN DIPOLE OF HALF-LENGTH h. 


#** REACTANCE IS IN THE LIMIT OF AN INFINITELY THIN FELEMENT- 


Table 4. Electrical Properties of Very Thin Monopole Elements on Groundplanes of Zero, 
Large, and Infinite Extent 


electrically-short and half-wave dipoles. The peak directivity of a 
quarterwave monopole is 1.88 dBi and 5.16 dBi for groundplanes of zero 
and infinite extent, respectively. The directivity on the horizon of 
a quarterwave monopole is 1.88 dBi and -0.86 dBi for groundplanes of 
zero and very large but finite extent, respectively. Slightly smaller 
directivities are obtained for electrically-short elements than for 


quarterwave elements. 


The radiation resistances obtained by different methods are 
compared in table 5 for a thin quarterwave element on a groundplane of 
radius 0 < ka < 8.5 wavenumbers. The suspected best available results 
are obtained by the integral equation method for 0 < ka ¢ 1.0, by the 
method of moments (N#1) for 1 < ka < 3.0, by the oblate spheroidal 
wave function method for 3.00 < ka ¢ 6.5, and by the method of moments 
(N=1) for 6.5 < ka < 14. The results obtained by the method of 
moments is in good agreement with the suspected best available results 


obtained by other methods. 


The numeric directive gain patterns of a quarterwave element on 
groundplanes of radii ka = 0, 3, 4, 5, J42, and ©0 wavenumbers are 
plotted in figure 8. These polar graphs of directive gain should not 
be confused with the polar graphs of radiated power density plotted in 
Ref. [5] for constant base current. It should be noted that the peak 
directivity and direction of peak directivity are not monotonic 


functions of the groundplane radius. 


The input impedance of a thin quarterwave element is plotted in 
figures 9 and 10 for groundplane radii 0 < ka < 14 wavenumbers. The 
input resistance varies between 19.4 and 46.1 ohms and asymptotically 
approaches 36.5 ohms for increasingly large groundplane radii. The 
input reactance varies between - 0 and +32.5 ohms and asymptotically 


approaches + 21.3 ohms for increasingly large groundplane radii. 
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Radiation Resistance (ohms) 


Norma lized 
Groundplane Radius ,| Integral Equation] Oblate Spheroidal Moment 
2x ald Wave Function Method Met hod 
0 *19.43 

0.25 *19.48 19.49 
0.50 *19.62 19.62 
0.75 *19, 86 19. 86 
1.00 *20.23 20.21 
1.25 20.76 *20.71 
1.50 21.51 *21.25 
1.75 22.59 *22.44 
2.00 24.15 *23..89 
2.25 26. 46 *25.99 
2.50 29.95 27.32 *29.02 
2.75 35.44 31.92 *33.24 
3.00 44.60 *37.48 38.62 
3.25 *43.01 44.12 
3.50 *46 06 47.57 
3.75 *45.55 47.35 
4.00 *42.67 44,43 
4.25 *39.23 40.58 
4.50 *36.23 37.13 
4.75 *34.00 34.46 
5.00 *32.57 32.68 
5.25 *31.93 31.70 
5.50 *32.13 31.53 
5.75 *33.23 32.26 
6.00 *35.23 34.04 
6.25 *37.85 36.94 
6.50 *40.33 40.56 
6.75 30.12 *43 53 
7.00 30.09 *44 20 
7.25 *42 .30 
7.50 *39.10 
7.75 *35.96 
8.00 *33.50 
8.25 *31.88 
8.50 *31.16 


*SUSPECTED BEST AVAILABLE RESULT 


Table 5. Radiation Resistance of a Thin Quarterwave Element at 
the Center of a Circular Groundplane of Radius 
O0<ka<¢8.5 Wavenumbers. (Sinusoidal current distribution 


assumed on clement.) 
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(b) 2na/A = 3 


(a) 277a/A=0 


6=0 
Z 
g =U 
2 
(c) 2na/k = 4 (d) 2na/k = 5 


Figure 8 


(f) 2na/r = o 


(e) 2na/\ =/42 


Elevation Directive Gain Patterns, for Any Azimuthal 
Direction, of a Quarterwave Element Mounted on a 
Groundplane of Radius a. (The patterns are polar graphs 


on the same linear scale.) 
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The directive gain on the horizon, peak directivity, and elevation 
angle of the peak directivity of a quarterwave element are plotted in 
figures 11 - 13 for groundplane radii 0 < ka < 14 wavenumbers. The 
directive gain on the horizon varies between +1.9 and -1.9 dBi and 
asymptotically approaches -0.86 dBi for increasingly large but finite 
groundplane radii. The peak directivity varies between +1.0 and 
+5.3 dBi and asymptotically approaches +5.2 dBi for increasingly large 
groundplane radii. The elevation angle of the peak directivity varies 
between 33 and 90 degrees and asymptotically approaches 90 degrees for 


increasingly large groundplane radii. 


The radiation resistances of thin elements of length h/A 20.25, 
0.1 and 0.025 for groundplane radii 0 < ka < 8 wavenumbers are 
compared in figure 14. The radiation resistance of each element is 
normalized to its radiation resistance for a groundplane of zero 
extent. The normalized radiation resistance, as a function of 
groundplane radius, is approximately independent of the element 
length. For an electrically-short, thin element (h/A < 0.1) whose 
length is small or comparable to the groundplane radius, the input 
reactance (not shown) is approximately independent of groundplane 


radius. 
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Figure 12. Peak Directivity of a Thin Quarterwave Element at the Center of a Circular 
Groundplane of Radius 0<¢ka <14 Wavenumbers 


14,0 


68 


90.0 


60.0 

ro] 

us 

o 

a 70.0 

ui 

at 

(C) 

4 

< 

z 60.0 

is} 

= 

< 

> 

rm 

a 

w 50.0 
40.0 
30.0 

0 
Figure 13. 


METHOD OF MOMENTS, N=1 
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NORMALIZED GROUNDPLANE RADIUS, 277a/A 


12.0 13.0 


Elevation Angle of Peak Directivity for a Thin Quarterwave Element at the Center 
of a Circular Groundplane of Radius 0 <ka< 14 Wavenumbers 
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SECTION 4 


MODELS IN WHICH THE CURRENT DISTRIBUTIONS ON THF MONOPOLF 
ELEMENT AND GROUNDPLANE ARE ROTH INITIALLY UNKNOWN 


4.1 Boundary Conditions 


In Section 3 the total field at an arbitrary field point could be 
expressed simply as the vector sum of an incident field and an induced 
field {see Eq. (3.1-12)] because the incident field was specified from 
the initially known current distribution on the monopole element. Tn 
this section, such a procedure is not possible because the current 
distributions on the element and groundplane are both initially 
unknown. Instead, in this section, the total field is determined by 
representing the unknown current distribution on either the element, 
groundplane, or both by a finite series of overlapping modes with 
unknown current amplitudes. The current amplitudes are determined hv 
matrix inversion subject to boundary conditions comprising current 


constraints and the excitation voltage across the coaxial] line feed. 


The antenna geometry is shown in figure 1 of Section 2.1. The 
current waveform is given by Eq. (2.1-2). The models which follow in 
this section are based on the current characterization and circuit 


idealization in figure 2(c) of Section 2.4. 


The coaxial line feed and excitation voltage, which is not 
explicitly shown in figure 1, is characterized in figure 2(c) by a 
surface magnetic current density (magnetic frill) Me which sits on top 
of a thick groundplane of radius b, where b, is the radfus of the 
outer conductor of the coaxial line feed. The magnetic frill Me is 


defined over the groundplane by Eq. (2.4-8). 
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Constraints on the various circuit currents of figure 2(c) are 
given by Eqs. (2.4-3) - (2.4-7). These constraints, together with the 
magnetic frill, constitute the boundary conditions on the current 


amplitudes. 
4.2 Method of Moments, 0 < ka ¢14 


When the element and groundplane current distribution are 
initially unknown, the current distributions may be determined by a 
sinusoidal-Galerkin moment method employed by Richmond‘), The 


antenna geometry is shown in figure 1 of Section 2. 


In the moment method, the element is divided into N equal segments 


(see figure 15) of length d” given by 
d“ = h/N, N is a positive integer (4.2-1) 


where h is the element length. The groundplane is divided into M 


concentric annular zones (see figure 16) of width d given by 
d = (a-b)/M > b, 7 b, Mis a positive integer (4.2-2) 


where a is the groundplane radius, b is the element radius, and b, is 


the radius of the outer conductor of the coaxial line feed. The : 
groundplane extends from the coaxial line inner conductor because of 
the equivalent circuit representation in figure 2(b) of the coaxial 
line excitation. The element current distribution I(z) and 
groundplane current distribution I(P) are the sum of the current 
distributions on each element segment and groundplane annular zone, 


respectively, and are given by 


N 
U2) = D> Lz) (4.2-3) 


n=1 
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Figure 15, 


P(r, 6.) 


Ina = MZ= 0) Q(0:77/2,.9') 


Element with N Equal Segments of Length d' 
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Figure 16. Groundplane with M Annular Zones of Width d 
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M 
1(P) = D> I2¢p) (4.2-4) 


n=] 


th 


where I (z) and TP) are the current distributions on the n- segment 


and oth annular zone, respectively. 


In the sinusoidal Galerkin method, the current distribution I (2) 


on the nth element segment Is approximated by 


I, sin[k(nd“-z)]} + I +l sin{k[z-(n-1)]d7} , 
sin(kd’ ) ? 
(n-1)d° < z < nd’; n#1,2,3, ....N 
(4.2-5) 


I (2) ~ 


where Li«1 = 1(z=0) and Ta = 0. 


The current distribution Th (P) on the at? 


by 


annular zone is approximated 


ae sin{k[md+b-p]} + V1 sin{k[0-(m-1) ]d-b} 


1 (P) = ——————sin(kdl) ; 
b + (m-1)d < P< btmd; m1, 2, 3, ...M 
(4.2-6) 
where lel = “Loi = -I(z20) and rey = 0. 


The current distribution on each segment and annular zone is 
therefore the sum of two overlapping dipole modes except for the 
segment and annular zone adjacent to the base of the element. The 
“base” mode may be thought of as a dipole mode comprising an element 
monopole mode (n=1) and a groundplane monopole mode (m=1) with 
terminals along the circumference of the base of the element. The 
element and groundplane dipole modes are functionally tied to each 


other by means of the base mode amplitude constraint Ii-l = -lhe1 al 
I(z=0). 
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The element and groundplane total current distributions are 
represented in Eqs. (4.2-3) - (4.2-6), as a series of N+M overlapping 
sinusoidal dipole modes with N+M unknown complex current amplitudes. 
The N+M amplitudes are determined in Richmond’s moment method by 
inversion of a (N+M) x (N+M) matrix subject to the boundary conditions 
discussed in Section 4.1. The constraint I-17 ag ee reduces the 
matrix size to (N+M-1) x (N+M-1). 


The numbers of subdivisions, N and M, are limited by the cost of 
computation time and by the precision of the computer. The accuracy 
of the solution can be decreased appreciably if either N or M is too 
small or too large. The method of moments converges to a solution 
when an increase or decrease of unity in the value of N or M does not 


appreciably alter the solution for input impedance. 


A method of moments computer program for a monopole element at the 
center of a circular groundplane in free space was obtained from 
Prof. Richmond of Ohio State University. The program computes the 
input impedance and the N+M complex current amplitudes on the element 
and groundplane for a voltage input V(O0) = 1 volt. The program, 
written in FORTRAN IV language and in single precision, was edited and 
converted by The MITRE Corporation to double precision for use on a DEC 
VAX 11/780 computer. A listing of the MITRE version, designated 
“RICHMD1”, is given in Appendix B2. 


Our experience with this program suggests the following 
constraints on the use of the program. At least double precision is 
required to give convergent results. Meaningful results were obtained 
for element radii b/A > 10°19 ana groundplane radii ka > 0.25. Asa 
rule of thumb, N#2-3 kh and M=2-3 ka give reasonably accurate and 
convergent results. The amount of central processing unit (CPU) time 


on the VAX 11/780 computer is approximately 2 minutes for N+M = 20, 6 
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minutes for N+M=30, 34 minutes for N+M=60, and 50 minutes for N4M=75. 
It is most likely that the CPU time could be considerably reduced if a 
more efficient method than the Crout method were employed for matrix 


inversion. 


The versatility of a sinusoidal - Galerkin method of moments is 
illustrated by the case N=l. A sinusoidal current distribution is 
imposed on the element by setting N=1 [see Eq. (4.2-5)]. 


The element and groundplane current distributions obtained from 
program RICHMD1 can be utilized to obtain the far-field radiation 


pattern in the following manner. 


The magnetic field intensity H(P) at a far-field point P(r, 6, 6) 
and for an e tut waveform is related to the magnetic vector potential 
ACP) py (26) 


— TaN — 

H(P) = ~j(1/M wu x A(P) (4.2-7) 
where a. is a unit vector in the radial direction and 7 is the 
free-space wave impedance. The magnetic vector potential A(P) at the 


far-field point P resulting from a current source point Q(x”, y”, 27) 


= Q(b, b, z”°) on the element or a source point Q(p”, 77/2, $6”) on the 


groundplane (see figure 4) is given bye 


AC) = (,/4™) . 5(Q) exp(-JkPQ) (1/PQ) dv 
Vv 


= [u, exp(-jkr)/47r] fs exp( jku_ - 09) dv (4.2-8) 
Vv 
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where 


Be = [rtnae O05 + xe tay Pee tee OG, n> (xrPay Zee 2y ll? 
a +0Q = x” sin@ cosd + y” sin® sind + z~ cos@ 


i permeability of free space = 47 x 1077 henry/m. 
J(Q) = source surface current density at an arbitrary source point Q 
on either the element or groundplane (amp/m) 
r = OP = radial distance between the origin 0(0,0,0) and the 
far-field point P(r,@ ,¢@) 
dv = differential surface area containing the current source 


points Q (m2) 


Substituting Eq. (4.2-3) into Eq. (4.2-8), the far-field magnetic 


field intensity a resulting from the monopole current distribution 


is given by 
L N 
Hy” = espe) sind f > 1_(2°) exp(jk2z°cos®@) dz” 
. 0 n=." 
j_exp(—jkr) 


N 
Um sind sin(kd7y LS I, exp( jnkd cos0)| 


n=l 


l-cos(kd“) cos(kd“cos®) - cos 8 sin(kd” )sin(kd“cos®) 
+ j cos(kd”) sin(kd” cos®@) - j cos @ sin(kd”~) cos(ka”cos@) | 


N 
a2 I, expl j(n-2) ka“cos0] | [cos(ka") cos(kd’cos @) 

=2 
+ cos 8 sin(kd”) sin(kd’cos 6) - 1 + j cos(kd”) sin(kd’cos 6) 
- j cos © sin(kd”) cos(kd’cos 0)| (4.2-9) 
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Substituting Eq. (4.2-4) into Eq. (4.2-8), the far field magnetic 
field intensity aay resulting from the groundplane current 


distribution is given by 


_ = a at M 
HA? i aLeycie) ff 3. x 1") exp[ jkP“ cos (6-6) } 


cos 8 cos(°-6) dd” dp” 

(4.2-10) 
Eq. (4.2-10) is not readily evaluated when Eq. (4.2-6) is substituted 
into it. An approximate simplified solution can be found when each 
annular zone of the groundplane is subdivided into X smaller annular 
zones of width A/P= d/X so that the current distribution in each 
subdivided zone is almost constant and approximately equal to its 
average current. The total number J of subdivided annular zones on 


the entire groundplane is given by 
J = XM (4.2-11) 


The average current I in the uth subdivided zone is given by 


u APtb 


ns Ae. 1a sln{k(md+b-p] + Ti+ sin{k[P- (m-1)d - b]} ie. 
u Ap (u-1)AP+ b sln(kd) 
[b+(m-1)d] <P<(b+md) (4.2-12) 


The current distribution I(f) on the groundplane is therefore given 
approximately by 


1(p) = )> i (4.2-13) 
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Substituting TP") for TC?) in Eq. (4.2-10), 


a a 


H(8) ~, J exp(~jkr) exp(= jkr) 
6 


exp(-jkr)cos 9 
47r sin® 


“ Yexp{ jkp*cos(°-0)]}cos® cos(¢°-)dd-do" 


J 
= Die") [sot sin@{(u-1)Ap +b]} - J, {ksin®[u ap+b])] 
ai (4.2-14) 


where Jy is the Bessel function of the first kind, of order zero. 


The total far-field magnetic field intensity aad is given by 
(total) (e) (g) 
H =H + 4.2-1 
6 6 * % arpo sree 
(e) (g) 


where Ny and Hg are given by Eqs. (4.2-9) and (4.2-14), respectively. 


The time-averaged radiated power density s(®), directive gain 
d(®), and radiation resistance R are found from Eqs. (2.2-4), 


(total) 


and (2.2-7) where Ny is given by Eq. (4.2-15) and I(z=0) = I 


n=1° 
Numerical evaluation of s(®), d(®), and R are implemented by means of 
MITRE computer program "RICHMD2" written in FORTRAN language for use 
on a DEC VAX 11/780 computer. A listing of program RICHMD2 is given 


in Appendix B2. 


A discussion of the results obtained by the method of moments, for 
a sinusoidal current distribution on the element (N=1), is given in 
Section 3.4. 


Computer printouts of the directive gain and far-field radiation 
pattern of the experimental monopole antennas discussed in Section 5 


were obtained by the method of moments and are given in Appendix A2. 
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The input impedance determined by the method of moments (N4, 
M=3ka) for a quarterwave element of radius b=10-°A and groundplane 
radii 6 < ka < 30 wavenumbers is compared in table 7 of Section 4.5 
with that obtained by the method of moments combined with the 
geometric theory of diffraction (GTD). For ka > 14, Richmond’s method 
of moments is inaccurate. In addition to the large CPU time required 
(greater than 30 minutes for ka > 14), convergent results are 
difficult to obtain for ka > 14. It is most likely that the 
usefulness of Richmond’s method of moments could be extended to 
significantly larger groundplane radii if a more efficient method of 
matrix inversion, different from the Crout method employed by 
Richmond, were utilized. When the element current distribution is not 
assumed, Richmond’s method of moments presently gives the best 


available results for groundplane radii ka < 14 wavenumbers. 


The lower limit of ka, for accurate results utilizing Richmond’s 
method of moments, has not been firmly established. For an element 
segment number N > 1, we have obtained useful results for ka as small 
as 0.25 wavenumbers but not for groundplanes of zero extent. For N=l, 
accurate results were obtained (subject to the constraint of a 
sinusoidal element current distribution) for groundplanes of zero 


extent. 


The radiation resistance of various diameter resonant (zero input 
reactance) elements is compared in figure 17 of Section 4.5 with the 
radiation resistance of an infinitely thin quarterwave element for 
groundplane radii 0 < ka < 8.5 wavenumbers. For ka —*0, the 
radiation resistance of the resonant elements is appreciably 
different from that of the infinitely thin quarterwave element because 
an infinitely thin element on a groundplane of zero extent is resonant 


for an element length equal to 0.5A and not 0.25, [see Eq. (32-38)]. 
The lengths and radiation resistances of various diameter resonant 
elements, for groundplanes of radii 0.25 < ka < 7 wavenumbers, are 


plotted in figures 18 and 19 of Section 4.5. 
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4.3 Method of Moments Combined with Geometric Theory of Diffraction, 
8.5 < ka <0 


The method of moments, when used to determine the current 
distributions on both the element and the groundplane, can require 
considerable computation time for large groundplane radii. For 
example, Richmond’s method of moments requires approximately one-half 
hour of CPU time on the VAX 11/780 computer for a groundplane radius 
ka=20 with N=4 element segments and M=60 groundplane annular zones 
(see Section 4.2). Although it may be possible to reduce the 
computation time by use of a more efficient program than the one 
employed by Richmond, it is of interest to find an alternative method 


for large groundplanes. 


Awadalla and Maclean(?)»(10) have reduced the computation time for 
large groundplanes by combining the method of moments with the 
geometric theory of diffraction (GTD). The element current 
distribution is determined by the method of moments and the effect on 
input impedance by groundplane edge diffraction is determined by GTD. 


The antenna geometry is shown in figure l of Section 2.1. The 
method of Awadalla and Maclean is described in [9] for determining the 
input impedance and in [10] for determining the radiation pattern. 


A fictitious magnetic edge current is defined and expressed in 
terms of a GTD diffraction coefficient to account for diffraction by 
tne edge of the groundplane. The method of GTD is valid only for 
sufficiently large groundplane radii. Unfortunately, this method, as 
applied by Awadalla and Maclean, does not determine the groundplane 


current distribution. 
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The groundplane current distribution is idealized by Awadalla and 
Maclean in determining the radiation pattern. The current 
distribution on the top of the groundplane is assumed to be that for 
an infinite groundplane but defined only over the finite extent of the 
groundplane. The net current at the edge of the groundplane is set 
equal to zero. This constraint satisfies the boundary condition given 
by Eq. (2.4-4) and requires, at the edge of the groundplane, that the 
current on the bottom of the groundplane be equal but in the opposite 
direction to that on the top of the groundplane. The current 
distribution on the bottom of the groundplane is then assumed to 
decrease linearly from the edge of the groundplane to zero at the base 


of the element. 


The method of moments employed by Awadalla and Maclean is not a 
sinusoidal-Galerkin method. Consequently, one cannot impose a 
sinusoidal current distribution on the element by setting the number 
of element segments N=l as a test case for purposes of comparison with 


other models. 


A program listing was obtained from Prof. Awadalla and edited at The 
MITRE Corporation. A listing of the edited program “AWADALLA” is 
given in Appendix A5. Program AWADALLA is written in FORTRAN language 
for use on a DEC VAX 11/780 computer. The CPU time for N=30 element 


segments is less than 10 seconds. 


Program AWADALLA was utilized to obtain the input impedance, 
directive gain, and far-field elevation pattern of a quarterwave 
element of radius b=107° on groundplane of radii 8 < ka < 50 
wavenumbers. The results for input impedance are given in table 6. 
These results are compared in table 7 of Section 4.5 with those 
obtained by the method of moments (N=4). Concerning input impedance, 
the method of moments combined with GTD gives useful results for 
ka > 6, accurate results for ka > 8, and the suspected best available 


results for 14 < ka <@. 
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TABLE 6. INPUT IMPEDANCE OF A QUARTERWAVE ELEMENT OF RADIUS b=1076 AT THE 
CENTER OF A CIRCULAR GROUNDPLANE OF RADIUS 8 < ka < 50 WAVENUMBERS 
(Method of moments combined with geometric theory of diffraction) 


INPUT RESISTANCE INPUT REACTANCE 
(OHMS) (OHMS) 


8 18.28 

9 24.45 
10 21.82 
11 18.67 
12 23.08 
13 22.51 
14 19.17 
15 22.18 
16 22.82 
17 19.70 
18 21.51 
19 22.89 
20 21.05 
22 22.84 
23 20.63 
24 20.73 
25 22.67 
26 21.06 
27 20.56 
28 22.46 
29 21.40 
30 20.47 
31 22.19 
32 21.68 
33 20.49 
34 21.93 
35 21.90 
36 20.56 
37 21.65 
38 22.06 
39 20.68 
40 21.42 
41 22.13 
42 20 .86 
43 21.21 
44 22.14 
45 21.05 
46 21.04 
47 22.10 
48 21.26 
49 20 .93 
50 21.99 
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A sample computer printout of directive gain and the far-field 
elevation pattern obtained by program AWADALLA for ka=20 is given in 
Appendix A5. Unfortunately, the idealizations, made by Awadalla and 
Maclean in characterizing the groundplane current distribution, yield 
unrealistic peak directivities and elevation patterns. For example, 
for ka=49, a peak directivity of 7.5 dBi was obtained. This result 
for a thin quarterwave element seems unlikely because a thin 
quarterwave element with a sinusoidal current distribution and mounted 
on a groundplane of infinite extent has a peak directivity of only 
5.2 dBi. 


The method of moments combined with GID via the definition of a 
fictitious magnetic edge current has also been reported by Thiele and 
Newhouse(!1) for computing input impedance and by Stutzman and 
Thiele28) for computing the far-field radiation pattern. In the 


(28) 


method of Stutzman and Thiele » the far-field pattern is determined 


without idealizing the groundplane current distribution. 
4.4 Method of Images, ka= co 


For the idealized case of a monopole element mounted on a 
groundplane of infinite extent and of infinite conductivity, a 
monopole antenna of length h may be modelled by the method of images 
as a dipole of total length 2h but with one-half the input impedance 
and double the peak numeric directivity of the dipole (see 
Section 3.8). The infinite groundplane prevents radiation into the 
hemisphere below the groundplane but generates fields in the upper 


hemisphere identical to those of a dipole. 


A detailed treatment of the fields and input impedance of a dipole 
is given in [13]. An excellent summary of the present 
state of the art of dipole theory, including plots of input impedance 
as a function of dipole length and dipole radius, is given in [14]. 
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4.5 Summary of Results 


In Section 4, the current distribution on the element is 
determined by the method of moments rather than being assumed as was 
the case in Section 3. The determination of the element current 


distribution is particularly important for element radii b > 107*). 


The essential difference, in the two models utilized in 
Section 4, is the treatment of the groundplane current distribution. 
In Richmond’s model, the groundplane current distribution is 
determined by the method of moments. Useful results are obtained for 
0 < ka < 14 wavenumbers. In the model of Awadalla and Maclean for 
large but finite groundplanes, the input impedance is accurately 
computed by the introduction of a fictitious magnetic edge current 
determined by the Geometric Theory of Diffraction (GTD). However, in 
that model the groundplane current distribution is not determined but 
instead is idealized -- causing unrealistic peak directivities and 


far-field radiation patterns. 


The input impedances, determined by the method of moments and the 
method of moments combined with GTD, are compared in table 7 for a 
quarterwave element of radius b=10 OA and groundplane radii 
6 < ka < 30 wavenumbers. The suspected best available results for 
input impedance are obtained by the method of moments for 0 < ka ¢ 8.5 
and by the method of moments combined with GTD for 8.5 ¢< ka < oo. The 
method of moments combined with GID is inaccurate for ka < 8 
(approximately 6% error in input reactance for ka=8 and 7% error in 
input resistance for ka=6). Richmond’s method of moments is 
inaccurate for ka > 14 (approximately 11% error in input reactance and 
5% error in input resistance for ka=15). Richmond’s method of moments 
would be useful for ka > 14 wavenumbers if a more efficient method 


than the Crout method were employed for matrix inversion. 
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TABLE 7. INPUT IMPEDANCE OF A QUARTERWAVE ELEMENT AT THE CENTER OF A CIRCULAR GROUNDPLANE OF RADIUS 
6 < ka < 30 WAVENUMBERS. (h/A= 0.25, b/A= 10 ') 


INPUT RESISTANCE (OHMS) INPUT REACTANCE (OHMS 


NORMALIZED (1) (2) C)=-(@) (2) (1)-(2)| 
GROUNDPLANE Moment Moment Method (1) x 100 Moment Method (1) x 100 
RADIUS, ka Method with GTD Method with GTD 


6 *35.2988 37.7258 6.88 *26.7927 26.6729 0.45 
7 *45.7499 42.6747 6.72 *20 .5664 20.2537 
8 *35.7335 35.9640 0.65 *17.1839 18.2823 
9 36.2886 *36 5496 0.72 24.7424 *24 4496 
10 42 .0987 *41 .4538 1.53 21.5758 *21.8225 
11 38.1401 *37 65375 1.58 18 .8073 *18.6728 
12 36.7590 *36 .2962 1.26 22.0824 *23 0841 
13 39.7689 *40.4928 1.82 22.6158 *22.5061 
14 39 .6986 *38 .4400 2.60 20.3449 *19.1676 
15 38.3512 *36 3630 5.18 20.0026 *22.1845 
20 38 .4250 *39 .3019 2.28 22.2706 *20.1706 
25 39.0594 *38.1746 2.26 21.5934 *22 6704 
38 .8703 *37 .8798 2.55 21.0052 *20.4663 


(1) RICHMOND, METHOD OF MOMENTS (N=4, M=3 ka) 
(2) AWADALLA AND MACLEAN, METHOD OF MOMENTS COMBINED WITH GEOMETRIC THEORY OF DIFFRACTION (N=30). 


*SUSPECTED BEST AVAILABLE RESULT 


The radiation resistance of various diameter resonant (zero input 
reactance) elements is compared in figure 17 with the radiation 
resistance of an infinitely thin quarterwave element for groundplane 
radii 0 < ka < 8.5 wavenumbers. The radiation resistance of the 
resonant elements is not too appreciably different from the 
quarterwave element for ka > 1 wavenumber. However as ka —+0, the 
radiation resistances of the resonant elements become increasingly 
different from that of the infinitely thin element because an 
infinitely thin element on a groundplane of zero extent is resonant 


for an element length equal to 0.5A and not 0.25, [see Eq. (3.2-38)]. 


The lengths and radiation resistance of various diameter resonant 
elements for groundplane of radii 0.25 < ka < 7 wavenumbers are 
plotted in figures 18 and 19. For these groundplane radii, the 
element resonant length hieg/ * varies from approximately 0.22 to 0.34 
wavelengths for element radii 107? < b/A< 1072 wavelengths. The 
resonant radiation resistance, for these groundplane and element 


radii, varies from approximately 21 to 65 ohms. 
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SECTION 5 


COMPARISON WITH EXPERIMENTAL RESULTS 


The input impedance and elevation patterns of several monopole 
antennas were measured, each at a different frequency within the 
frequency band 30-254 MHz, on the MITRE Corporation VHF antenna range. 
Each antenna consisted of a tubular element of radius b=0.25 in. 
mounted at the center of a groundplane of radius a=4 ft. and fed by a 
50 ohm RG-214 coaxial cable with a type N panel connector and a 50 ohm 
tapered adapter to the element. The length h of each element was 
approximately a quarterwave. The exact length of each element length 
was chosen to be resonant (input reactance = zero ohms) for a 
groundplane of infinite extent. On a groundplane of finite extent, 
the input reactance is expected to asymptotically approach zero as 


ka —> 00. 


The VHF antenna range is located on the roof of MITRE E-Building. 
The transmitting and receiving antennas are at a height 27 ft. above 
the roof and are separated by approximately 40 ft. A conducting fence 
16 ft. high and 48 ft. wide is located on the roof midway between the 
transmitting and receiving antennas to minimize multipath reflections 
from the roof. Lossy ferrite toroids (Ferronics, Type B material, 
0.540 in. I. D. x 0.870 in 0. D. x 0.25 in.) are spaced along the 
receiving and transmitting cables to minimize currents on the outside 
of the cables. Outside rf interference is reduced to at least 40 dB 


below the desired signal by the use of narrowband rf filters. 
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The monopole test antenna was operated in the receiving mode to 
obtain elevation patterns and on the same platform for input impedance 
measurements. The monopole groundplane was supported by a 10 ft. 


wooden vertical mast mounted on an antenna pedestal. 


The input impedance was measured with a Hewlett-Packard HP-8754A 
Network Analyzer with a computerized printout. The measurement test 
set-up is shown in figure 20. The predicted input impedance of each 
experimental monopole antenna was determined by Richmond’s method of 
moments utilizing program RICHMD1. The predicted and measured input 
impedances are compared in Table 8. The measured input resistance 
differs from predicted values by approximately 1-10% over a range of 
normalized groundplane radii = 0.77 - 6.5 wavenumbers. The measured 
input reactance differs by approximately 2-12 ohms from the predicted 
values. Since some of the predicted values are near resonance, a 


percentage comparison is not made for input reactance. 


Measured elevation patterns of most of the test monopole antennas 
specified in Table 8 are compared in figures 21 - 29 with theoretical 
patterns predicted by the method of moments (see Appendix A2 for the 
computer printouts). Allowing for some multipath distortion by the 
VHF test range, the received patterns are in good agreement with the 
predicted pattern. It should be noted the measured pattern for ka = 
0.766 is not appreciably different from that predicted for the 
monopole element itself. The effect of the groundplane on the shape 


of the pattern is not readily apparent until ka 2 2.0. 


110 


OF TWICE THE LENGTH 
OF TEST CABLE 


HP8748A 


8-PARAMETER 
TEST SET-UP 


HP6328A 


FREQUENCY 
COUNTER 


HP8764A 


NETWORK 
ANALY ZER 


HP8601A 


STORAGE 
NORMALIZER 


HP2673A 


PRINTER 


HP9 836 


TEST CABLE 


COMPUTER 


LAN GUAGE 
AND DATA 
CONTROL 
OIsCcs 


DEVICE 
UNDER 
TEST 


Figure 20. Input Impedance Test Set-up 


lll 


<«@— REFERENCE PLANE EXTENSION CABLE 


BASIC LANGUAGE 
BASIC LANGUAGE EXTENSIONS 
MITRE PROGRAM ’NTWKANAL’ 


rau 


Feequeacy 
(M2) 


tines 


0.9 
1393.43] 
36.0 
[327.86 | 
43.0 
[274.48] 
54.0 
218.57] 
62.4 
(189.15] 
75.0 
157,37] 
86.0 
137.24] 
89.7 
[131.58] 
97.5 
121.05] 

117.0 
107.88} 
136.5 
[86.47] 
156.0 
[75.66] 
175.5 
[67.25] 
195.0 
(60.53) 
214.5 
[55.02] 
| 234.0 
(50.44 | 
253.5 
[46.56 


ee Xa 


Element 
Leagth , 
tin.) 


| normalized 
h 


i Groundp lane 
f Kadius, a4 
\ 


$ 
i 


Normalized T 
blemene ' 
Radius , F 
(b/d) x 10” 


Normalized 

Element 

Length , 
hid 


0.2382 
0.2379 
0.2374 
0.2366 
0.2366 
0.2363 
0.2355 
0.2352 
0.2346 
0.2344 
0.2335 
0.2330 
0.2328 


0.2335 


Input Resistance, Rj (chms) _|, Input Reactance*, Xig (ohms) | 

Predicted By Measured T Percentage 1, Predicted By Measured Dif ference 

Moment Method By MITRE**! Difference |, Moment Method By MITRE*™| (3) - (4) 
kL (2) Q)-@) x100 (3) (4) (chas ) 


* The »lement lengths were selected to be resonant when mounted on a yroundplane of infinite extent [see R.W.P. King, Theory of Linear Anteanas (Harvard 


University Press, 1956) 


tig. 30.7a). 


On -roundplanes of :inlte extent, the input reactance 


** MITRE measurements were made on !~-7-85 and 1-10-85 utilizing computerized printout of network analyzer. 


Table 8. 


Input Impedance of 0.50 Diameter Monopole Elements on an 8 ft. Diameter 
Groundplane. 


Is »xpected to asymptotically approach zero as ka —* oo. 


(Comparison of moment method predictions with MITRE measurements.) 


Figure 2]. 


z way od} 


+ RELATIVE POWER ON 


Measured and Theoretical Elevation Patterns, 
(h/A =0.2396, b/A =6.35 x 1074, ka = 0.766 


30 MHz 


FATTERN NO Of 
_.. [PROMCT 9600 
FNGRE 


Remcan, WS 


Ax 


OATE g@-v7-bai 


‘wOMOPOLE 
000 wits 
te ashwe ves 
ort matoe 
a Oe FT. MULTIOdTH 8CaR: 
+O FY, OTAINWELE Bow 


DATE @-a7 


PATTERN NO 8 
-{PROJEC! §ecoe 


ENGRS 


anxs M4 MOnoro.e 


Oe asAs8 ote 


40 71 manag 


eosats 


FY MUL TIPAT OC 


ovr 


acne 


erainweL 


i \ 
eo @¢ & 7 
‘AW® 3NO B3MO4 SAIL 


oa 
4 


' 
Lv 

! j- ' , 1 
| 


i 
9 7] 7 aan 


. I + - 
AVM 3NO -Y3MOd BAUYIIY! 


ee 


ae ae, Si 
W380 = 
7 


0.919) 


on Patterns, 36 MHz 


» Ka 


ti 
=7.63 x ioe 


Measured and Theoretical Elev 


(h/A =0.2396, b/A 


Figure 22. 


; 1 if : H 
: : i ’ 7 H 2 | : | 
+9 > ow SF Nw a2 + © GF a + w&© & 
cacy INO €3VGd BAU via. , : 7 , an us 1 
‘ ‘ H a - a 

[ihtet = IO Gi niet Aon | 

yt oat aaa eS 
: : f \ t 
: i : 3 , ! \ ; i 


zg 
tT 
wy 
arn 
YW OV 
cm 
wo 
(a 
wr 
ea) 
@ Il 
Au 

0 
cx 
io) 
ed a 
ast 
ca) 
Po 
ver 
et 
<a tad 
ams 
O <tr 
YU e 
‘do 
PS anal 
M il 
4 
On 
o™ 
oO 
a 
IN 
Go 
oo 

N 
ma 
vo 
ot 
a 
oO 
y 
= 
st 
Oy 


13.21 x 1074, ka = 1.593) 


0.2379, b/A = 


Measured and Theoretical Elevation Patterns, 62.4 MHz 


(h/A 


Figure 25. 


Measured and Theoretical Elevation Patterns, 75 MHz 


(h/A =0.2374, b/A = 15.89 x 1074, ka = 1.915) 


Figure 26. 


Oa BAe a eas a evan aean | 


TT, 6D AVM BND NAMOd GALY BE, ; 
Tee ca ae ay en ee ae | 


: ciel mtn Ey ah a oar 


eae aie ie 
: 
a 
{ 
| 
| 
= 
| 
| 
a 


a ines 
Nae PP: aries 


‘ . 


1 . i 
i : : 


1} 
|: 
i 4 
{ 
H 


86 MHz 


18.22 x 1074, ka = 2.197) 


and Theoretical Elevation Patterns, 


(h/A =0.2366, b/A 


Figure 27. Measured 


0.2355, b/A -24.78 x 1074, ka = 3.000) 


Measured and Theoretical Elevation Patterns, 117 MHz 


(h/d 


Figure 2g. 


i Tea eee 
Fe canst Sa eat ane ORG eee ee oe | 


| ea, ak 


“i a ae 
" et Sed ier oe! 


= 4.000) 


etical Elevation Patterns, 156 MHz 
or *, ka 


ed and Theor 
(h/A =0.2346, b/A =33.04 x 1 


29. Measur 


Figure 


10. 


ll. 


REFERENCES 


J. Bardeen, "The Diffraction of a Circularly Symmetrical 
Electromagnetic Wave by a Co-axial Circular Disc of Infinite 
Conductivity," Physical Review, Vol. 36, pp. 1482-1488, November 1930. 


J.H. Richmond, "Monopole Antenna on Circular Disk," Technical Report 
711639-1, Electro Science Laboratory, Ohio State University, July 
1979. See also IEEE Transactions on Antennas and Propagation, Vol. 
AP-32, No. 12, pp. 1282-1287, December 1984. 


A. Leitner and R.D. Spence, "The Oblate Spheroidal Wave Functions," 
Franklin Institute Journal, Vol. 249, No. 4, pp. 299-321, April 1950. 


A. Leitner, "Effect of a Circular Groundplane on Antenna Radiation," 
Mathematics Research Group, Research Report No. EM-19, New York 
University, Washington Square College, April 1950. 


A. Leitner and R.D. Spence, "Effect of a Circular Groundplane on 


Antenna Radiation," Journal of Applied Physics, Vol. 21, pp. 
1001-1006, October 1950. 


C.L. Tang, "On the Radiation Pattern of a Base-Driven Antenna Over a 
Circular Conducting Screen,” Technical Report No. 301, Cruft 
Laboratory, Harvard University, May 1959. 


J.E. Storer, "The Impedance of an Antenna Over a Large Circular 
Screen," Journal of Applied Physics, Vol. 22, pp. 1058-1066, August 
1951. See also Technical Report No. 119, Cruft Laboratory, Harvard 
Unviersity, November 1950. 


J.E. Storer, "The Radiation Pattern of an Antenna Over a Circular 
Ground Screen," Technical Report No. 126, Cruft Laboratory, Harvard 
University, June 1951. 


K.H. Awadalla and T.S.M. Maclean, "Input Impedance of a Monopole 
Antenna at the Center of a Finite Groundplane,"” IEEE Transactions on 
Antennas and Propagation, Vol. AP-26, No. 2, pp. 244-248, March 1978. 


K.H. Awadalla and T.S.M. Maclean, "Monopole Antenna at the Center of a 
Circular Groundplane: Input Impedance and Radiation Pattern,” IEEE 


Transactions on Antennas and Propagation, Vol. AP-27, No. 2, pp. 
151-153, March 1979. 


G.A. Thiele and T.H. Newhouse, "A Hybrid Technique for Combining 
Moment Methods with the Geometrical Theory of Diffraction," IEEE 
Transactions on Antennas and Propagation, Vol. AP-23, No. 1, pp. 
62-69, January 1975. In a telephone conversation with M. Weiner on 
July 28, 1983, G. Thiele reported that his computer program for a 
monopole on a circular groundplane, which he wrote while at Ohio 
State, was lost in moving to the University of Dayton. 


123 


12. 


13. 


14. 


15. 


16. 


17. 


18. 


19. 


20. 


21. 


22. 


23. 


24. 


25. 


26. 


27. 


28. 


REFERENCES (Concluded) 


R,F. Harrington, Time-Harmonic Electromagnetic Fields, New York: 
McGraw-Hill, 1961, pp. 106-112. Please note that the signs of the 
radial field Ey on p. 111 and the magnetic frill Mg on p. 112 differ 
from those in Eq. (2.4-9) and (2,4-8), respectively, because in 
Reference [12] the voltage at the outer conductor is positive with 
respect to that at the inner conductor. 


R.W.P. King, The Theory of Linear Antennas, Cambridge, MA: Harvard 
University Press, 1956, p. 87. 


R.S. Elliot, Antenna Theory and Design, Englewood Cliffs, New Jersey: 
Prentice-Hall, 1981, pp. 290-291. 


C.A. Balanis, Antenna Theory: Analysis and Design, New York: Harper 
and Row, 1982, p. 315. 


Op. cit. 15., p. 101, Eq. (4-2). 


L. Brillouin, "Origin of Radiation Resistance," Radioelectricité, Vol. 
3, pp. 147-152, 1922. 


J.A. Stratton, Electromagnetic Theory, New York: McGraw-Hill, 1941, 
pp. 454-457. 


M. Abramowitz and I.A. Stegun, Handbook of Mathematical Functions, New 
York: Dover Publications, 1972, pp. 231 and 244. 


Op. cit. 15, p. 292. 

W. Magnus and F. Oberhettigner, Formulas and Theorems for the 
Functions of Mathematical Physics, New York: Chelsea Publishing Co., 
1954, p. 34. 


I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Series and 


Products, New York: Academic Press, 1965. 


J.B. Thomas, An Introduction to Statistical Communication Theory, New 
York: John Wiley & Sons, Inc., 1969, Appendix B: The Dirac Delta 
Funct ion. 

Op. cit. 19, p. 763. 

Op. cit. 15, pp. 118-126, 133-135, 285-295. 

Op. cit. 14, p. 29. 

Op. cit. 15, pp. 502-522, 529-531. 


W.J. Stutzman and G.A. Thiele, Antenna Theory and Design, New York: 
John Wiley & Sons, 1981, pp. 489-495, 500-507. 


124 


